The following article is Open access

Multiphase Gas Interactions on Subarcsec Scales in the Shocked Intergalactic Medium of Stephan's Quintet with JWST and ALMA

, , , , , , , , , , , , , , and

Published 2023 July 6 © 2023. The Author(s). Published by the American Astronomical Society.
, , Citation P. N. Appleton et al 2023 ApJ 951 104 DOI 10.3847/1538-4357/accc2a

Download Article PDF
DownloadArticle ePub

You need an eReader or compatible software to experience the benefits of the ePub3 file format.

0004-637X/951/2/104

Abstract

We combine James Webb Space Telescope (JWST) and Hubble Space Telescope imaging with Atacama Large Millimeter Array CO(2–1) spectroscopy to study the highly turbulent multiphase intergalactic medium (IGM) in Stephan's Quintet on 25–150 pc scales. Previous Spitzer observations revealed luminous H2 line cooling across a 45 kpc-long filament, created by a giant shock wave, following the collision with an intruder galaxy, NGC 7318b. We demonstrate that the Mid-Infrared Instrument/F1000W/F770W filters are dominated by 0–0 S(3) H2 and a combination of polycyclic aromatic hydrocarbon and 0–0 S(5) H2 emission. These observations reveal the dissipation of kinetic energy as massive clouds experience collisions, interactions, and likely destruction/recycling within different phases of the IGM. In 1 kpc-scaled structure, warm H2 was seen to form a triangular-shaped head and tail of compressed and stripped gas behind a narrow shell of cold H2. In another region, two cold molecular clumps with very different velocities are connected by an arrow-shaped stream of warm, probably shocked, H2 suggesting a cloud–cloud collision is occurring. In both regions, a high warm-to-cold molecular gas fraction indicates that the cold clouds are being disrupted and converted into warm gas. We also map gas associated with an apparently forming dwarf galaxy. We suggest that the primary mechanism for exciting strong mid-IR H2 lines throughout Stephan's Quintet is through a fog of warm gas created by the shattering of denser cold molecular clouds and mixing/recycling in the post-shocked gas. A full picture of the diverse kinematics and excitation of the warm H2 will require future JWST mid-IR spectroscopy. The current observations reveal the rich variety of ways that different gas phases can interact with one another.

Export citation and abstract BibTeX RIS

Original content from this work may be used under the terms of the Creative Commons Attribution 4.0 licence. Any further distribution of this work must maintain attribution to the author(s) and the title of the work, journal citation and DOI.

1. Introduction

Since the discovery of a filament of radio continuum in the intergalactic medium (IGM) of the Stephan's Quintet (Allen & Hartsuiker 1972), this compact galaxy group has been studied at many wavelengths to try to better understand the remarkable nature of its multiphase IGM. As suspected in early studies (Moles et al. 1997; Sulentic et al. 2001; Xu et al. 1999), the overall picture that has emerged is that one of the galaxies, NGC 7318b, is colliding with the diffuse IGM of the main group at a very high velocity, creating a giant (45 kpc) filament of shocked gas seen from X-rays (Trinchieri et al. 2005; O'Sullivan et al. 2009), in the UV (Xu et al. 2005), and in ionized gas emissions (Xu et al. 2003; Iglesias-Páramo et al. 2012; Konstantopoulos et al. 2014; Duarte Puertas et al. 2019, 2021). Studies of the stellar populations along the giant shocked structure also suggest that star clusters are beginning to form there (Gallagher et al. 2001; Fedotov et al. 2011), although at a low rate, perhaps because some of the forming clusters lie inside bubbles of highly excited ionized gas (Konstantopoulos et al. 2014).

The physical conditions within the gas along the main shock front are far from simple and require more detailed study. Mid-IR spectroscopy of the main filament has showed that the entire filament is radiating strongly in pure rotational lines of molecular hydrogen (see Figure 1; Appleton et al. 2006; Cluver et al. 2010). An explanation of the excitation of such strong, L(H2) > 1041 erg s−1, molecular hydrogen emission from a region experiencing fast shocks from the intruder's high velocity was presented by Guillard et al. (2009). The model assumes that the main intruder shock propagates into a clumpy pre-shock medium, heating low-density regions to X-ray temperatures, but causing mildly overdense regions to collapse to form H2 on grains on timescales shorter than the dynamical timescale of the collision. These molecular clouds experience low-velocity molecular shocks as they continue to dissipate mechanical energy from their surroundings. The interaction of these warm H2 clouds with their surroundings is a key area that we will investigate with the James Webb Space Telescope (JWST) and Atacama Large Millimeter Array (ALMA) in this paper.

Figure 1.

Figure 1. Stephan's Quintet: a gray-scale representation of the NIRCam F150W ERO image of Stephan's Quintet with overlays of the Spitzer IRS image of the 0–0 S(1) H2 line (blue contours) and the half-power primary beam size of the ALMA CO(2–1) observations (red circles) described in Section 2.2. The faintest features seen in the F150W image are at a surface brightness level of 0.7–1 MJy sr−1. Contours of H2 emission are in units of 0.3, 0.53, 0.75, 0.98, 1.2, 1.43, 1.65, 1.88, and 2.1 MJy sr−1 from Cluver et al. (2010).

Standard image High-resolution image

A more detailed analysis of the H2 excitation properties of the warm H2 in Stephan's Quintet was discussed in Appleton et al. (2017), where it was shown that in order to heat so much H2 (109 M) at low temperatures (150–400 K), a mix of low-velocity magnetic C-shocks (∼5–10 km s−1) and faster (15–25 km s−1) J-shocks were needed.

Very broad spectral line widths have been observed in the Stephan's Quintet filament. Guillard et al. (2012) used the IRAM 30 m telescope to detect CO emission from several regions in the main filament and the bridge, and revealed line profiles of several hundred kilometers per second in at least three separate kinematic components. In addition to gas detected at the systemic velocity of the intruder NGC 7318b (5774 km s−1), and that of the remaining group members (6600 km s−1), broad lines of CO emission were also detected at intermediate velocities, suggesting that molecular gas has formed out of diffuse shocked gas. Recent UV spectroscopic observations with the Hubble Space Telescope's (HST) Cosmic Origins Spectrograph (COS) targeted five regions along the main filament and in the bridge, and detected extremely broad Lyα emission (FWHM exceeding 1500 km s−1) over small sampled regions of ∼1 kpc scale (Guillard et al. 2022). The lines were even broader in some cases than those seen in the [C ii] line (Appleton et al. 2013) using Herschel, suggesting some resonant scattering of UV photons. The detection of broad-line, powerful Lyα emission, broad [C ii] neutral gas, and broad CO emission strongly supports the picture of a highly turbulent medium containing many gas phases resulting from a turbulent cascade of energy from the large to the small scales.

Although the Early Release Observations (ERO) made by JWST of Stephan's Quintet cover the entire inner group (including NGC 7319, NGC 7318a/b, NGC 7317, and the foreground galaxy NGC 7320), this paper will concentrate on the main shocked filament and bridge that were observed by Spitzer in spectral mapping mode (Cluver et al. 2010; Appleton et al. 2017) using the InfraRed Spectrograph (IRS; Houck et al. 2004). As shown in Figure 1, this included the main north–south H2 filament between NGC 7318b/a and NGC 7319, as well as the H2 bridge connecting NGC 7319 to the main filament. It is known that much of the gas, including H i, ionized gas, and molecular gas, lies outside the main body of the galaxies. Although the previous IRS observations of the mid-IR lines covered much of this gas, it did not cover all. The Long-Low IRS coverage included all of the area shown in Figure 1, but IRS Short-Low coverage was restricted to only the the main H2 filament and part of the H2 bridge (see Section 3.1).

Data from the JWST and ALMA allow us to directly compare, for the first time at comparable resolution, the distribution of warm 0–0 S(3) H2 (as we will demonstrate via the F1000W MIRI band) to cooler molecular hydrogen mapped through the low-J CO(2–1) line, as well as ionized gas emission (Hα with WFC3 HST, and Paα with NIRCam F200W). Very hot gas, which forms an important IGM component in Stephan's Quintet, is detected in X-rays throughout, and especially in the main shock. Its distribution cannot be compared in detail to our current observation, because even the deepest Chandra images (O'Sullivan et al. 2009) do not detect enough photons at arcsecond scales to allow meaningful comparison.

These observations will help us gain an understanding of how turbulent energy, driven mainly by the intruder galaxy, is dissipated from large driving scales to smaller dissipation scales, and how this affects the cooling of the gas through different gas phases and temperature regimes; and, in particular, how can we explain so much energy flowing out through low-velocity shocks needed to produce the high radiant line luminosity from warm molecular hydrogen discovered by Spitzer. Although we do not expect to probe down to the smallest dissipation scales, we will show that we already see differences in the distribution of different thermal phases in the IGM at ∼150 pc scales probed by ALMA and JWST.

The current paper will not try to address all aspects of the extended emission detected by JWST, but will primarily concentrate on two main objectives. First, we compare the JWST MIRI images with the Spitzer IRS spectra and spectral images to identify the dominant features present in the three MIRI filters (F770W, F1000W, and F1500W). Second, we compare the distribution of warm molecular gas, cold molecular gas, and ionized gas in three contrasting regions within the main IGM between NGC 7319 and NGC 7318b/a. The three regions we have chosen to emphasize in this paper were observed with ALMA in CO(2–1). The primary beams (fields of view) of the three ALMA pointings and their associated field numbers are also identified in Figure 1. 18 We will discuss the regions in more detail in Section 4. Although we touch on the sparsity of recent star formation in two out of three of the studied fields, we leave a full discussion of the star formation properties and cluster formation to a later paper.

In Section 2, we will describe observations made with JWST, HST, ALMA, and the Combined Array for Research in Millimeter-wave Astronomy (CARMA). Given that we do not yet have spectroscopy of the IGM in Stephan's Quintet, we present, in Section 3, a discussion of what each of the MIRI and NIRCam images is likely to contain by comparing the JWST images with Spitzer spectral maps of especially warm H2 and polycyclic aromatic hydrocarbon (PAH) bands. In Section 4, we present the results of our zoomed-in study of three major emission complexes in the IGM observed at high resolution by ALMA in the CO(2–1) line. This allows us to synthesize from JWST, ALMA, and HST the best possible information we have about the nature of the warm and cold molecular gas, the ionized gas, and the formation of star clusters, all at subarcsec resolution. Section 4 also includes a discussion of the relative fraction of warm and cold molecular gas obtained for the regions using both the JWST and ALMA data. In Section 5, we discuss the results in the context of our understanding of turbulence and shocks, which we believe largely dominate the gas dynamics of this multiphase IGM. In Section 6, we present our conclusions.

We assume in this paper a distance to Stephan's Quintet of 94 Mpc (for H0 = 70 km s−1 Mpc−1 and an assumed group heliocentric systemic velocity of 6600 km s−1) for consistency with previous work (e.g., Appleton et al. 2017; Guillard et al. 2022). At this distance 1'' corresponds to a linear scale of 456 pc. The overall scale of the giant north–south intergalactic filament in Stephan's Quintet is 47 kpc.

2. Observations

2.1. JWST Images

The first description of the NIRCam and MIRI ERO observations, and the choice of Stephan's Quintet as a target (PID 2732), is presented in Pontoppidan et al. (2022). The NIRCam images were taken in six filters (F090W/F277W, F150W/F356W, F200W/F444W), with an integration time of approximately 20 minutes for each band. The FULLBOX five-point "5TIGHT" dither pattern was used, resulting in a large rectangle mosaic of three pairs of dithered tiles, covering 6.3 × 7.3 arcmin2. The MIRI image covers a much smaller field of view than the NIRCam mosaic, i.e., the central galaxies, NGC 7318a/b, NGC 7319 (a Seyfert 2 galaxy), and NGC 7320 (foreground galaxy), using four tiles in three MIRI bands, F770W, F1000W, and F1500W. The integration times were about 22 minutes per filter. The large cycling eight-point dither pattern was used to maximize the spatial coverage of a single tile.

We worked with the level 2b images retrieved from MAST. All the HST and JWST data used in this paper can be found in MAST: doi:10.17909/dfsd-8n65.

2.2. ALMA Observations

Observations of CO(2–1) in the three fields of Stephan's Quintet were made with the ALMA 12 m array in a mosaic observing mode on 2015 July 21 (ID: 2015.1.000241.S; PI: P. Guillard). The total integration time for the mosaic was 35 minutes. The spectral setup consisted of four spectral windows of 1.875 GHz with 3.9 MHz channels. One of the spectral windows was centered on the redshifted CO(2–1) (νrest = 230.538 GHz), while the remaining three spectral windows covered line-free continuum emission. The standard ALMA calibration plan included bandpass, phase, and flux calibration.

The data were calibrated with the scriptForPI.py calibration script that was included with the archival data, using the ALMA calibration pipeline version r32491 that is included with the Common Astronomy Software Applications (CASA) version 4.3.1 (CASA Team et al. 2022). After calibration, the continuum emission was subtracted in the (u,v) domain by fitting a straight line to the line-free channels. The line data were subsequently imaged using the CLEAN algorithm to produce data cubes with 20 km s−1 channels and a resolution of 0farcs36 × 0farcs21 (PA ∼ 10°). The rms noise in these data cubes is 0.5 mJy beam−1 chan−1.

A primary beam correction was applied to produce accurate flux measurements. The half-power beamwidth of the primary beam is 22farcs8 (10.4 kpc). Total intensity maps of the CO(2–1) emission were made by integrating the signal across the channels where line emission was detected.

2.3. Combined Array for Research in Millimeter-wave Astronomy CO Observations

Observations of Stephan's Quintet were made with CARMA on 2010 August 1 (program c0593; PI: K. Alatalo) using 15 antennae in the CO(1–0) line. The original data were processed with MIRIAD (Sault et al. 1995). Data reduction and imaging were done in identical fashion to Alatalo et al. (2013). The original data were cleaned (Högbom 1974), which resulted in a restored synthesized beam of 4.1 × 3.3 arcsec2 (1.9 × 1.5 kpc). Velocities from extracted spectra have been converted to heliocentric velocities assuming a shift of −11 km s−1 from kinematic local standard of rest to heliocentric. Velocities are all quoted using the optical definition of velocity.

2.4. HST Archival Images

We dearchived calibrated images from the HST MAST archive for filters F665W and F336W with the WFC3/UV instrument. For F665N, the image encompasses the Hα line and includes data based on five dithered observations, each taken with integration times of 5200 s each in 2009 July and August, and obtained as part of a WFC3 ERO campaign (SM4/ERO). For F336W, the observations were obtained as part of proposal ID 12301 (PI: S. Gallagher) on 2011 October 27, with a total integration time of 14,474 s. The resolution of the WFC3/UV observations is comparable with that achieved by NIRCam (0farcs05–0farcs06). For F665W and F336W, the expected FWHM of a point source is 0farcs07, and 0farcs075, respectively, corresponding to ∼30 pc at D = 94 Mpc.

2.5. WCS Alignment of the JWST Images

The observations were obtained from the JWST and HST MAST archive, and most images required small adjustments to the WCS coordinates for proper alignment. We initially used Gaia Data Release 3 (DR3) stars to align one image (the HST WFC3 F665N) to the DR3 system using the STScI software package WCSTweak. The NIRCam images were also similarly aligned. However, for MIRI the same method did not produce good results, probably because of the smaller number of stars used for alignment. We used a refinement scheme to solve this problem as we need good alignment between all the images. This was especially important for one of the target fields, which required continuum subtraction of the F1000W and F770W images with F1500W. First, we created small subimages of the F1000W, F770W, and F1500W images around the Field 6 area. The subimage contained at least three Gaia DR3 stars close to the target of interest. These subimages were then aligned more carefully by making small pixel-level shifts in each image to bring the subimages into close alignment with each other and the positions of the Gaia stars. We applied the same method to Fields 5 and 4 (see Section 4). The final results produced local MIRI images which aligned to better than 1 MIRI pixel (0farcs3) over the small scale of the extracted regions. The results were validated where the CO ALMA images aligned with clearly related features.

Following Papadopoulos et al. (2008), the uncertainty in the absolute astrometry of the ALMA data is δ θbas = (δ B · Δ k )/ B ≈ (δ ϕbas/2π)〈Θbeam〉. The phase error δ ϕbas ∼ (2π/λ)(δ B · δ k ), with ∣δ B ∣ ∼ 0.1 mm being the assumed typical error in the baseline length, ∣δ k ∣ = 4.26° the distance to the phase calibrator J2216+3518, λ ∼ 1.33 mm the wavelength of the redshifted CO(2–1) emission, and Θbeam ∼ 0farcs36 the major axis of the synthesized beam. This results in δ θbas ∼ 0farcs02, which means that the uncertainty in the absolute astrometry of the ALMA data is small compared to that of the HST and JWST data. 19

3. Expectation of the Emission Observed through the MIRI and NIRCam Filters in Three Study Regions

This section will explore our expectations for what emission processes are likely to contribute to the MIRI images based on previous Spitzer spectral mapping results. This paper will concentrate mainly on three regions which lie near the centers of the ALMA CO(2–1) observing fields. The regions are typical of the diversity of mid-IR spectral properties encountered in the IGM of Stephan's Quintet (SQ) in general, and are the focus of this paper. They are highlighted in Figure 2 as boxes (dotted lines) superimposed on one of the publicly released ERO images of the inner part of the group. Field 6, which we will show is dominated by 10 μm H2 emission (green color in this figure), was targeted by HST COS (Guillard et al. 2022) and shows extremely broad Lyα emission, broad [C ii]158 μm emission and contains some of the the warmest H2 observed by Spitzer (e.g., Cluver et al. 2010). It is representative of a shock-heated region in the main north–south filament. Field 5 is another H2-dominated region which lies outside the main north–south H2 filament, and observations at other wavelengths (e.g., single-dish CO observations by Guillard et al. 2012, and Herschel [C ii] observations by Appleton et al. 2013) show that this gas is also highly turbulent. It forms part of the bridge of H2 emission seen by Spitzer. Field 4, by contrast, lies in a well-studied extragalactic star-forming region, SQ-A (sometimes referred to as the northern star-forming region) and has been studied extensively at optical wavelengths (Gallagher et al. 2001; Xu et al. 2003). Its star-forming nature can be inferred from its bluish-white appearance in the false-color map of Figure 2 due to a combination of weaker H2, strong PAH 7.7 μm emission, and dust continuum (see Cluver et al. 2010; Appleton et al. 2017).

Figure 2.

Figure 2. Stephan's Quintet: a false-color representation of the MIRI F770W, F1000W, and F1500W JWST image of the inner Stephan's Quintet, showing the three regions highlighted in the discussion that are detected near the primary beam centers in the ALMA CO(2–1) emission-line observations. Field 6 is representative of the center of the main filament, where previous observations with Spitzer have shown that the molecular hydrogen has the highest temperature. Field 5 lies in the bridge region between NGC 7319 and the main shock identified in Figure 1. Finally, in Field 4, a bright region in the SQ-A star-forming region is highlighted (see text).

Standard image High-resolution image

3.1. Comparison of JWST MIRI Images with Spitzer IRS Spectral Mapping Results

Under conditions found in the disks of normal galaxies, strong 7.7 μm and 8.6 μm PAH features are expected to dominate the MIRI F770W filter for low-redshift galaxies (e.g., Smith et al. 2007). In Stephan's Quintet (Appleton et al. 2006, 2017; Guillard et al. 2009; Cluver et al. 2010) much of the mid-IR spectrum of the IGM is dominated by strong pure rotational H2 lines and emission from [Si ii]λ34.8 μm. Unlike normal star formation regions, the majority of the main IGM filament in Stephan's Quintet exhibits very weak 7.7 μm PAH complex emission and weak nebular lines like [S iii]λ18.7, 33.5 μm (Cluver et al. 2010). This is the result of low star formation rates in the main filament (measured to be between 0.05 and 0.08 M yr−1; Cluver et al. 2010), and emission lines more typical of fast shocks than H ii regions (e.g., Konstantopoulos et al. 2014). Examples of these H2-dominated spectra have been presented in Figures 8, 13, and 14 of Cluver et al. (2010) and Figures 15 and 16 of Appleton et al. (2017). The Spitzer IRS spectra showed that only a small minority of places along the main filament (including parts of the SQ-A region) are spectra found with line flux ratios more typical of star-forming regions.

In Figures 3(a)–(c), we compare the filter transmission bandpass for the MIRI imaging bands (F770W, F1000W, and F1500W) directly with the Spitzer IRS spectra taken within Fields 4, 5, and 6 (Figure 1). The MIRI filter F770W is likely to contain emission primarily from the 7.7 μm PAH complex in regions where star formation dominates, but in other regions, especially those with strong H2 emission, the 0–0 S(5) can also potentially contribute, as in Figure 3(a). The F1000W filter, on the other hand, almost exclusively contains line emission from the S(3) mid-IR pure rotational H2 line (see also Figures 15(e) and (h) of Appleton et al. 2017). Based on the presented spectra, we expect the MIRI F1500W filter to be dominated by faint dust continuum from the general IGM, as well as from warm dust from embedded star clusters and sometimes the extended wing of the broad 17 μm PAH complex (Smith et al. 2007); see center panel of Figure 3.

Figure 3.

Figure 3. Spitzer IRS spectra of three 5.5 × 5.5 arcsec2 (2.5 × 2.5 kpc2) regions centered on the three fields studied at high resolution in CO(2–1) emission with ALMA and at 7.7 and 10 μm with JWST MIRI imaging (see Figure 5). The MIRI broadband filter transmissions are superimposed to emphasize that the F1000W MIRI filter detects almost exclusively emission from the S(3) line and that the 7.7 μm band captures a combination of the weak 7.7 μm PAH complex and emission from S(5) H2. The 15 μm MIRI bandpass detects mainly faint dust continuum and some (weak) [Ne iii] emission. ALMA Field 5 was not covered by the Short-Low IRS module. In regions of strong star formation the F1500W filter can also be contaminated with the wings of the 17 μm PAH complex (Smith et al. 2007).

Standard image High-resolution image

A clean separation between the mainly PAH-dominated emission in the F770W MIRI filter and the 0–0 S(3) H2-dominated F1000W band can be seen in many regions along the main IGM filament and in the bridge region in Figure 4. Here, we overlay the contour maps obtained by isolating a PAH band 20 from the IRS spectral map from Spitzer in Figure 4(a) (white contours) with a two-color image of the JWST MIRI F1000W (green) and the MIRI F770W (orange). In Figure 4(b), we overlay the 0–0 S(3) H2 contours (red) over the same image. The overlays show that the white contours of PAH emission from Spitzer follow quite closely the brighter emission in the F770W-dominated regions, whereas the red contours of H2 follow closely the green emission originating mainly from the F1000W filter. This figure demonstrates how the 10 μm MIRI filter, as we expected from the individual IRS spectra, detects preferentially H2 emission, while the F770W filter detects those regions dominated by PAH emission. As shown in Cluver et al. (2010), the PAH emission is relatively faint in the main shock, but this figure emphasizes those regions with stronger PAH emission. This near isolation of the H2 emission in many regions of the IGM is a result of the relatively narrow width of the F1000W filter and the lack of other lines or bands that can seriously contaminate it for the redshift of Stephan's Quintet.

Figure 4.

Figure 4. Stephan's Quintet: a false-color representation of the MIRI F770W and F1000W JWST images with contours of rest-frame (a) 11.3 μm PAH and (b) 9.66 μm emission derived from spectral cubes obtained from the Spitzer IRS Short-Low mapping. The white and red solid lines show the extent of the IRS spectral mapping (see Cluver et al. 2010). The figure demonstrates how the S(3) line isolated in the Spitzer cube follows closely the 10 μm MIRI emission (green), and the PAH features follow the MIRI 7.7 μm image that contains contributions from the 7.7 μm PAH and H2 emission (see also Figure 3 of this paper and Appleton et al. 2017).

Standard image High-resolution image

3.2. Dust Continuum

Observations using the F1500W MIRI filter are expected to detect mainly warm dust emission, similar to that mapped by Spitzer at 24 μm with MIPS (see Guillard et al. 2010). However, in the SQ-A region (Figure 3(b)), the dust continuum is also contaminated by [Ne iii] emission. As noted above, the very strong 0–0 S(1) line is, fortuitously, redshifted out of the F1500W band.

Before discussing results for each of the zoomed-in regions in the next section, we mention the removal of faint dust emission from F1000W and F770W at the center of Field 6 (near the center of the main IGM filament in Stephan's Quintet). In this field, faint dust continuum in the F1500W filter was seen in the vicinity of the structures we are interested in (see Section 4). We therefore removed the faint dust continuum from both the F1000W and F770W images, using the Spitzer IRS spectrum of the region as a guide. First, we created small subimages of the field in all three bands centered on the main structure of interest and corrected them for small offsets in the WCS, as discussed in Section 2. After removing a median sky background from each subimage, we then subtracted the F1500W image from both the F770W and F1000W images, scaling the continuum by a factor of 0.9 to account for a slight rise in continuum seen between 7.7, 10, and 15 μm in the IRS spectrum of the region. This allowed us to better define the distribution of the H2 and PAH features. The effect was quite small (the dust emission was quite weak compared with the emission seen in the 7.7 and 10 μm bands) on the resulting fluxes of the F1000W and F770W band. Without spectroscopy of the region on the same spatial scale as the observed features, this method is necessarily imperfect, and where we estimate the fluxes for the features observed in Field 6 in the next section we include larger uncertainties to take this into account. For the subregion near the center of Field 5, no obvious dust structure was seen in the F1500W band near the structures of interest, and so we did not attempt to continuum subtract the H2 and PAH-dominated emission from the F1000W and F770W images. For Field 4, centered in the SQ-A star-forming region, our attempts at subtracting a scaled version of the dust continuum led to a severe oversubtraction of the emission in the F770W and F1000W bands for most reasonable dust scale factors. This is likely due to a combination of factors. First at that position (see Figure 3(b)), we have shown that the 15 μm band shows emission from not only dust continuum but also relatively strong [Ne iii], as well as broad PAH-band emission from the 17 μm "plateau" complex (see Smith et al. 2007). These additional contaminants make it difficult to isolate the warm dust component in the 15 μm image. Therefore, for this star-formation-dominated field, we conclude that the F1000W and F770W MIRI filters are hopelessly contaminated by dust emission in a way that cannot easily be corrected without high-resolution spectroscopic data from the MIRI Medium Resolution Spectrometer (MRS). This significantly limits our discussion of the molecular hydrogen properties in Field 4. For that region, we limit our discussion of the warm H2 properties determined from the Spitzer IRS observations.

3.3. Near-IR Images

We also present images in this paper obtained through the NIRCam F200W filter. F150W provides a near-IR continuum band relatively clear of emission lines at the redshift of Stephan's Quintet (0.022 < z < 0.025), with the well-known near-IR [Fe ii] lines falling just blueward or redward of the filter transmission. The bandpass of the F200W filter is sensitive to emission from the hydrogen recombination line Paα and starlight. Diffuse emission from this line becomes obvious when we compare some of these images with those obtained in the F150W NIRCam filter, and by comparison with the HST Hα images.

JWST/MIRI, which is able to detect warm H2 through the pure rotational transitions, can significantly improve our knowledge of how energy is dissipated in the IGM down to the 0farcs3 scales of giant molecular cloud (GMC) complexes (∼140 pc). NIRCam, with its ability to probe down to scales of 0farcs05–0farcs06 (23–27 pc at 1.5 and 2 μm, respectively) can probe ionized gas and stellar associations and clusters at scales at which sporadic star formation along the main filament is observed.

4. Results

4.1. Analysis of the Three Regions

In this section, we explore a comparison of archival HST and JWST MIRI data with the ALMA CO(3–2).

4.1.1. Field 6: Within the Main Shock-dominated North/South Filament

In Field 6 (Figure 5(a)), we show the integrated CO(2–1) emission from cold molecular gas superimposed on the warm H2 emission. The warm H2 gas (assumed to dominate the F1000W filter from the 0–0 S(3) H2 line) is distributed in an elongated structure extending over 4''–5'' (∼2 kpc) with a triangular-shaped hot spot to the southwest and a series of fainter clumps extending to the east. The overall distribution of warm H2 has the appearance of a head-tail structure. The CO emission (white contours) displays a clumpy, narrow shell-like structure associated with the warm H2 head, and a scattering of CO clumps along the tail extending over the same overall extent but not correlating in detail with the hot spots in the warm H2.

Figure 5.

Figure 5. A zoom in on the center of the ALMA Field 6 (see Figure 1) comparing (white) contours of CO(2–1) emission to (a) continuum-subtracted pure warm 0–0 S(3) H2, (b) continuum-subtracted warm 0–0 S(5) H2 plus 7.7 μm PAH complex emission, (c) Hα emission from F665N WFC3 HST, (d) mix of IR (red compact) stellar sources and Paα emission in the NIRCam f200w filter, (e) blue star clusters from F336W WFC3/UV HST (Gallagher et al. 2001). In (f), we show a composite RGB color representation of the F1000W (mainly warm H2), ALMA CO(2–1) emission (cool H2), and the F200W (Paα plus star clusters) images. For clarity, this subfigure does not overlay the CO contours. The MIRI identifications of the contributing lines/PAH bands are clear from the Spitzer IRS spectra of a region 5.5 × 5.5 arcsec2 (2.5 × 2.5 kpc2) in Figure 3(a), which spatially covers this entire region. All images were carefully aligned using Gaia DR3 stars close to the region shown, leading to relative uncertainties in astrometry of ∼0farcs15. The F770W and F1000W JWST images were smoothed to the same resolution as the F1500W image before subtracting the carefully aligned dust continuum. The white scale bar is 3'' (1.4 kpc) in length. A more detailed description of the CO emission and its kinematics is given in Figure A1.

Standard image High-resolution image

Based on the published IRS spectra, Figure 5(b) is interpreted as a mix of H2 (0–0 S(5)) and weak PAH emission. The dominant triangular-shaped head structure seen in the 10 μm image is also strongly represented here. A second major clump of 7.7 μm emission is seen to the northwest, which is weak at 10 μm. It is plausible that this feature is PAH dominated, whereas the emission from the compact head (which is seen also at 10 μm) may be dominated by warmer than normal H2, which would emit strongly in the 0–0 S(5) line. Without mid-IR spectroscopy, this cannot be verified.

This distribution of ionized gas can be inferred from the narrowband HST F665N filter images (dominated by Hα emission) in Figure 5(c) and the F200W NIRCam image (dominated by Paα emission and near-IR starlight) in Figure 5(d).

Both images show significant extended (presumed ionized gas) emission associated with the head of the warm H2 including a protrusion, or "spike," labeled in Figure 5(c). The feature is also seen in CO emission. Extended (ionized) gas also follows more closely the distribution of the warm H2 in the tail than it does the CO emission. This is well demonstrated in the composite image Figure 5(f), which shows an RGB representation of the 10 μm (warm H2), CO emission (cold H2), and F200W (ionized gas and stars) images. Gemini optical spectroscopy (Konstantopoulos et al. 2014), which covered this region at a spatial resolution of 1'' (458 pc), strongly suggested that the ionized gas is shock-excited, which would explain why the ionized gas follows the warm H2 emission in the head and the clumpy tail, and why it is less correlated with the CO emission.

The spatial resolution of the F200W image is far superior to the HST image, and it is clear that there are several compact sources on scales of <0farcs05–0farcs6 (23–27 pc) embedded in the ionized gas component. These sources are likely stellar associations or unresolved super star clusters (SSCs; Gallagher et al. 2001). They are generally rather sparsely distributed and poorly correlated with the warm or cool molecular hydrogen.

One exception is a group of compact sources and extended F200W emission associated with the southern tip of the northwest clump of CO emission identified in Figure 5(c). These sources may be evidence of recent star formation associated with the head-tail structure. Evidence of this comes from the bright u-band sources associated with the same region shown in Figure 5(e). Although a full analysis of the NIRCam colors of the point sources in Stephan's Quintet will be discussed in a separate paper, we know these blue regions have the optical colors of young clusters with ages estimated to be between 3 and 5 Myr based on previous HST UBVmVI photometry (Fedotov 2014; see Table B1 for more details). The existence of young star formation near the northwest CO clump may also explain why that region exhibits PAH emission, since PAHs are believed to be excited by UV radiation from young stellar associations.

The kinematics of the diffuse ionized gas centered on the head from optical spectroscopy shows extremely broad line widths (>800–1200 km s−1) in several key emission lines, e.g., [O iii]5007, [O i]6300, and Hα (Konstantopoulos et al. 2014; Guillard et al. 2022). The broad line widths have been attributed to a combination of turbulence and large-scale motions associated with gas along the line of sight. Even broader wings in the Lyα profile were also seen at this position. We will argue in Section 5 that the warm H2 emission detected by JWST at this position may be responsible for UV scattering of Lyα emission within the turbulent regions and, if so, should exhibit broad line widths.

While we cannot yet measure the kinematics of the warm gas, the CO observations provide kinematic information about the cold molecular phase. Figures 6(a) and (b) show the integrated CO surface density (moment 0) map and the mean velocity field (moment 1) map of the CO emission, respectively. The systemic velocity of the CO emission falls between the velocity of the intruder galaxy NGC 7318b (V(int) = 5770 km s−1) and that of the velocity barycenter of the main group (V(group) = 6600 km s−1). This is consistent with gas which has been decelerated in the collision of the intruder with the group-wide gas. However, the line width of the entire region in CO is less that 120 km s−1, as shown in the integrated ALMA CO(2–1) spectrum of the entire region and the CARMA CO(1–0) spectrum shown in Figure 6(c). This width is significantly smaller than that seen in the ionized gas, implying that the cold molecular gas does not take part in the turbulent motions seen in the ionized gas phase. More extensive kinematic information for the CO emission is presented in Appendix A, where we show individual spectra extracted from many regions (Figure A1), along with tabulated properties of the CO emission on the scale of the ALMA beam (Table A1). The tip of the head structure has the highest radial velocity, and the clumps show a trend to lower velocities as one moves east into the tail, reaching the lowest velocity in one of the most easterly clumps. We will argue later that this suggests material is being stripped from the head into the tail through ram pressure stripping with a hot medium.

Figure 6.

Figure 6. ALMA Field 6 region (labeled in Figure 2): (a) CO(2–1) integrated emission, (b) the intensity-weighted mean velocity map of the same region, and (c) the ALMA CO(2–1) and CARMA CO(1–0) spectrum of the region 5'' × 5'' (2.5 × 2.5 kpc2) centered on the same region. Effective synthesized beam shape (0.36 × 0.21 arcsec2 = 164 × 96 pc2 FWHM) is shown graphically in the left-hand corner of each figure. Arrows indicate the radial velocity of the intruder V(int) and the barycentric velocity of the main group V(gr).

Standard image High-resolution image

The velocity dispersion 21 of the individual clumps around the structure show a mix of unusually broad line widths (up to 100 km s−1 FWHM) and narrow lines (<40 km s−1). In Table A1, we show that the region near the "spike" feature seen in the ionized gas is unusually broad (FWHM =102 km s−1). Away from the head, several regions in the tail of the structure show line widths ranging from 60 to 90 km s−1 (FWHM; see Figure A1) to values lower than 40 km s−1. Given the masses of the individual molecular clumps of a ∼few ×106 M (Table A1), it is unlikely that the clumps with high-velocity dispersion are gravitationally bound. 22 The broader lines are consistent with turbulent motions on the subarcsec (∼100–150 pc) scale in the colder gas component, while others have much narrower lines.

One such region occurs in the CO emission clump closest to the bright blue star clusters described earlier. Here, the CO emission is essentially unresolved (40 km s−1 or less), perhaps suggesting turbulent motions are calmer there. Table A1 presents information about the line fluxes and estimated H2 masses for all the regions measured in detail with ALMA.

4.1.2. Field 5: A Shocked Structure in the Bridge Region

Field 5 lies outside the main molecular filament in Stephan's Quintet, and is closer to the large face-on galaxy NGC 7319. It forms part of an apparent bridge of H2 emission discussed by Cluver et al. (2010). The JWST images (e.g., Figure 4) show that it is the most easterly of a series of irregular shock-dominated (appearing green in that image) clouds that run approximately east–west across the field and may not be physically connected. In contrast to Field 6, which contained a single coherent CO structure, Field 5 contains three separate bright structures visible in the CO maps (labeled in Figure 7(a)). These consist of a clumpy CO ring about 1'' across (458 pc), a compact elongated core 1farcs5 to the southeast of the ring, and a broken linear filament of gas running north–south 0farcs5 further east of the compact core. The direction of the elongation of the core points toward the center of the ring. Both Figures 7(a) and (b) show that the compact CO core lies near the brightest emission at the center of the arrow-shaped 10 μm (warm H2) structure. As the arrow-shaped emission narrows down toward the northwest, it connects to the northern part of the CO emission from the ring. Figure 7(c) shows diffuse Hα emission from the compact CO core but little emission from the ring. The higher spatial resolution of the NIRCam F200W image shows that the corresponding Paα in the core is also elongated Figure 7(d). Faint F200W emission is also seen in the brightest CO ring clump. The U-band image, Figure 7(e), shows no obvious emission near any of the CO molecular structures. Although not shown here, the F150W NIRCam image shows point-like sources in the northern part of the CO ring, suggesting that the F200W image is revealing a mix of H ii regions and compact sources within the ring. The composite image, showing the three gas phases in one figure, emphasizes the warm H2 connection between the ring and the core (Figure 7(f)).

Figure 7.

Figure 7. A zoom in on the center of the ALMA Field 5 (see Figure 1) comparing (white) contours of CO(2–1) emission to (a) pure warm 0–0 S(3) H2, (b) warm 0–0 S(5) H2 plus 7.7 μm PAH complex emission, (c) Hα emission from F665N WFC3 HST, (d) mix of IR (red compact) stellar sources and Paα emission in the NIRCam F200W filter, (e) blue star clusters from F336W WFC3/UV HST. In (f), we show a composite RGB color representation of the F1000W (mainly warm H2), ALMA CO(2–1) emission (cool H2), and the F200W (Paα plus star clusters) images. For clarity, this subfigure does not overlay the CO contours. The MIRI identifications of the contributing lines/PAH bands are likely from the Spitzer IRS spectra of a region 5.5 × 5.5 arcsec2 in Figure 3, which in this case only includes the longer IRS Long-Low wavelengths. All images were carefully aligned using Gaia DR3 stars close to the region shown, leading to relative uncertainties in astrometry of ∼0farcs15. The F770W and F1000W JWST images were smoothed to the same resolution as the F1500W image before subtracting the dust continuum. The white scale bar is 3'' (1.4 kpc) in length.

Standard image High-resolution image

A string of clumpy sources at 2 μm is clearly seen approximately 1'' to the east and north of the ring (F200W image). It is not clear if these sources have any real connection to the CO structures, as some could be background galaxies. One of the clumps has a faint U-band association (see Appendix B and Table A3), suggesting recent star formation. Near-IR spectroscopy may help to determine if they are physically associated with Stephan's Quintet by providing evidence of H ii-region-like emission lines at the Quintet's redshift, such as Paα.

The kinematics of CO gas in the compact core is strikingly different from that of the ring and the linear structure. Figures 8(a) and (b) show the integrated and radial-velocity map of the CO(2–1) emission, and spectra of the three different CO-emitting structures are presented in Figure 8(c). The velocities within the partial ring increase clockwise around the ring from 6631 km s−1 in the north, reaching the highest value in the south and east quadrant of 6688 km s−1. In contrast, the compact core (which lies near the peak in the warm H2 emission) has a much lower (negative 250 km s−1) radial velocity than the ring (∼6400 km s−1). The core also shows a velocity shear along its length of about 30 km s−1 over 0farcs5 (200 pc), as shown in Appendix A. Finally, the third linear structure to the east of the compact core shares velocities similar to that of the ring, with velocities ranging from 6620 km s−1 in the south to 6680 km s−1 in the north.

Figure 8.

Figure 8. ALMA Field 5 region: (a) the CO(2–1) integrated emission at a resolution of 0.36 × 0.21 arcsec2, and (b) the intensity-weighted mean velocity maps of the same region, (c) the ALMA CO(2–1) spectra of the three different extracted regions shown as red polygons in (a). Note the very different systemic velocity of the compact structure B. In (d), we overplot the CARMA CO(1–0) over a 5 × 5 arcsec2 aperture (blue line), compared with the velocity of the ring C. The effective ALMA beam shapes (FWHM ellipses) are shown graphically in the bottom-left corner of the images. The black vertical line on the spectra show the barycentric systemic velocity of the V(group).

Standard image High-resolution image

An in-depth view of the kinematics of Field 5 in given in Appendix A, including individual spectra, Figure A2, and a table of CO properties, Table A2. The spectra show that individual CO clumps in all of the three CO structures exhibit broad CO line widths at the scale of the ALMA beam (∼150 pc). The compact core itself has a velocity dispersion of 145 km s−1 (some of this may be due to the gradient seen along the core), and all components of the ring have line widths which lie in the range 60–120 km s−1 (FWHM). The linear structure shows narrower line profiles in the range 48–87 km s−1. Appendix A also shows that, although the majority of the gas in the ring has much higher radial velocities than the core structure, one small segment of the ring (at the point nearest the core) has a velocity similar to the core. This is supported by the CARMA spectrum, Figure 8(d), which shows a blueward component associated with the average spectrum of the ring. This strengthens the idea that the two CO structures are somehow related despite their discrepant radial velocities. We will discuss in the next section the possibility that the warm bridge connecting the CO ring and core may be splashed warm molecular gas caused by the collision of a dense molecular cloud moving at high (blueward) velocity with respect to more diffuse gas associated with material in the group.

Some of the CO emission in Field 5 may be associated with a larger filament and other CO clumps, as seen in a larger-scale (10''–15'' = 4.5–7 kpc) image, Figure 9(a), where we present the ALMA CO emission map smoothed (to 1 × 1 arcsec2) to bring out fainter features. This reveals a very faint filament of CO gas and several fainter southern CO concentrations. Both the filament and the southern CO structures have faint warm H2 counterparts. Filaments like this one, in this case extending over 6''–10'' in scale (3–5 kpc), seem to be common in the MIRI maps of SQ (e.g., Figure 2). Figure 9(b) shows again the velocity field of the ring, compact core, and linear structure, this time in relation to the larger-scale structure in the region. This further emphasizes how the compact core has such a different velocity from all the other structures in the region. The velocity of the filament is poorly determined in CO because it is so faint, but seems to have a velocity similar to that of the average velocity of the group. Mid-IR spectroscopy would be needed to provide more information about the filament since it is well detected in the MIRI 10 μm image (presumed H2 emission).

Figure 9.

Figure 9. Extended CO(2–1) emission for ALMA Field 5 over a larger area than Figure 7. In (a), the CO emission (white contours) has been smoothed to a resolution of 1 × 1 arcsec2 (450 × 450 pc2) and superimposed on the warm 0–0 S(3) H2 detected within the F1000W filter. In (b), the smoothed version of the velocity field is shown for the same region. This emphasizes the unusually low radial velocity for the compact core compared with the other surrounding regions. The long filament connecting the brighter CO structures is seen in the warm H2. The filament's radial velocity is shown schematically because it is weak, but has the same velocity as the majority of the other structures except the compact core. The white scale bar in (a) is 3'' (1.4 kpc) in length.

Standard image High-resolution image

4.1.3. Field 4: In the Northern Star-forming Region, SQ-A

Finally, we present the distribution of gas phases and star clusters in the Field 4 SQ-A region. The region has been discussed previously in terms of a possible tidal dwarf galaxy (Moles et al. 1997; Mendes de Oliveira et al. 2001), although this is disputed by Xu et al. (1999), who point out that the kinematics do not show simple rotation. Figure 10 shows a closer correspondence between the various gas phases, compared with Fields 5 and 6. Panels (a) and (b) of Figure 10 show a close correspondence between the CO distribution, the 10 μm H2 image, and the 7.7 μm PAH-dominated image. In some of the panels, we have smoothed the CO distribution to an effective resolution of 1'' (458 pc). The close correspondence between these three interstellar medium (ISM) indicators is not unexpected since previous IRS spectroscopy of the region showed that the total H2-to-PAH ratios in that region are more typical of normal photo dissociation region (PDR) regions associated with H ii regions. Also, ground-based spectroscopy shows that the optical emission-line ratios from the ionized gas are consistent with H ii-region-like spectra (Konstantopoulos et al. 2014). Extended patchy ionized gas is seen in panels (c) and (d) of Figure 10. A sharp boundary in the Hα distribution near the point where the CO emission starts to rise on the eastern edge, suggesting significant extinction of Hα light due to a dust lane associated with the denser molecular gas. This is supported by the fact that the same "edge" is not seen in the F200W NIRCam filter (Figure 10(d)), where the Paα is likely to be much less extinguished that the Hα. Konstantopoulos et al. (2014) measured high optical extinction in this region based on ground-based optical spectra. The emission seen in the F200W filter, Figure 10(d), shows many slightly extended clump sources scattered across an elongated "disk." The overall scale of this possible starburst dwarf system is 5 × 1.5 arcsec2 (2.3 kpc × 700 pc). The U-band image supports the idea of strong extinction in the region of the peak CO emission. In Table B1, we also present UBVmVI photometry of some of the brighter well-defined sources (Fedotov 2014). This study shows that the brighter clumps have ages based on the optical photometry ranging from 6 to 36 × 107 yr, much longer than that of the blue clumps we investigated in Field 6. This may suggest a more mature stellar population, although because of higher extinction here NIRCam photometry would provide a more definitive result.

Figure 10.

Figure 10. CO(2–1) emission (white contours) at the center of ALMA Field 4 (the SQ-A northern star-forming region; see Figure 1) at two different spatial resolutions. The CO emission in (a), (b), (c), (e), and (f) has been smoothed to a resolution of 1 × 1 arcsec2 (456 × 456 pc2), whereas in (d) the full resolution (0.36 × 0.21 arcsec2 (164 × 96 pc)) CO emission is shown. CO contours are superimposed on (a) the band containing the warm 0–0 S(3) H2 line, (b) the band dominated by 7.7 μm PAH complex emission for these northern star-forming regions, (c) Hα emission from F665N WFC3 HST, (d) mix of IR (red compact) stellar sources and Paα emission in the NIRCam f200w filter, (e) the blue star clusters from F336W WFC3/UV HST. In (f), we show a composite RGB color representation of the F1000W (mainly warm H2), ALMA CO(2–1) emission (cool H2), and the F200W (Paα plus star clusters) images. For clarity, this subfigure does not overlay the CO contours. Unlike Figures 5 and 7, because of strong contamination of the F1500W continuum band by [Ne iii] emission (see Figure 3(b)) no continuum subtraction was made for the MIRI images for (a) and (b). All images were carefully aligned using Gaia DR3 stars close to the region shown, leading to relative uncertainties in astrometry of ∼0farcs15.

Standard image High-resolution image

The structure of the CO is that of two interleaved half-shells (very well seen in the smoothed version of the CO contours, which are shown in all the panels of Figure 10 except (b) and (d)). As with the previous descriptions, Figure 10(f) is an RGB composite of the F1000W (red, mainly H2 dominated), the smoothed ALMA CO(2–1) (green, cool H2), and the F200W (red, Paα plus starlight) images to show the spatial relationship between the different gas phases and star clusters.

The kinematics of the system led Mendes de Oliveira et al. (2001) to identify the region as a potential tidal dwarf because of its clear rotation using Fabry–Pérot observations with the Canada–France–Hawaii Telescope. They derived a systemic velocity of 6680 km s−1, an asymmetric rotation curve, and a high-velocity dispersion of 115 km s−1, which may have been the result of a blend of two different optical emission-line components.

Our CO spectra do not show obvious rotation. We show the integrated CO(2–1) emission map (moment 0) and velocity field map (moment 1) for Field 4 in Figure 11. Although the overall kinematics of the CO gas in Field 4 shows lower systemic velocity on the eastern side (6700 km s−1) and a higher value to the west (6740 km s−1), the velocity range is narrow (FWHM = 47 km s−1), centered on the barycentric velocity of the main SQ group. The overall systemic velocity is in reasonable agreement with the optical emission-line velocities, although there is little real hint of rotation. This may argue against the formation of the object within a tidal stream, where the rotation is expected as a remnant of tidal shear. As we show in Table A3 and Figure A3, the velocity dispersion of the clumps at the highest resolved scale is narrow, which is very different from the more turbulent emission in Fields 5 and 6.

Figure 11.

Figure 11. ALMA Field 4 in the SQ-A region: (a) CO(2–1) integrated emission associated with the brighter lower-velocity (Comp. 1) in the center of ALMA Field 4 associated with the star clusters shown in Figure 10, (b) the intensity-weighted mean velocity map of the same region, and (c) the ALMA CO(2–1) and CARMA CO(1–0) spectrum of the region 5'' × 5'' centered on the same field. Note that both ALMA and CARMA show a higher-velocity CO component (Comp. 2) associated with the region which is faint and diffuse in the ALMA observations. (a) and (b) only show the structure for Comp 1. Effective synthesized beam shape (0.36 ×0.21 arcsec2 FWHM) are shown graphically in the left-hand corner of each figure.

Standard image High-resolution image

Although not shown in Figure 10, there is a second very faint ALMA CO velocity component (6900 km s−1), which is seen at a higher amplitude in CO(1–0) observed with the larger beam of CARMA data. It is very clear as a double-horned profile in the CARMA spectrum of Figure 11(c). This emission must be quite diffuse and is only just detected with ALMA. A similar high-velocity CO component was seen associated with SQ-A by Guillard et al. (2009) with the 30 m IRAM telescope. There is a hint of this higher-velocity component in the kinematic maps obtained in Hα by Duarte Puertas et al. (2019), where the ionized emission is quite localized to SQ-A. So, like Field 5, we find two CO components associated with a major CO concentration in SQ with quite discrepant velocities. In this case, the two components are roughly coincident and so its difficult to know whether the warm H2 gas, which dominates the 10 μm JWST image, originates in two clouds in physical contact or are simply chance associations. Only mid-IR spectroscopy will allow us to fully explore this possibility.

4.2. Mid-IR Line Luminosities, Warm and Cold H2 Masses

In Table 1, we present measured and derived properties of the warm H2 emission derived from both the Spitzer IRS data of Appleton et al. (2017) and from measured fluxes from the F1000W image for Fields 5 and 6, where we are confident that the emission is primarily from the 0–0 S(3) H2 line. Unlike the tables in the Appendix A, which are concerned with illustrating the detailed kinematics and cold gas properties at the scale of the CO clumps, here we are concerned with the masses of the warm H2 and the cold molecular mass integrated over the same larger areas, presented in Table 2. For Field 4, where we were unable to properly correct for underlying dust emission across the map, we only quote the properties of the H2 from the IRS data on a large scale for comparison with Fields 5 and 6.

Table 1. Measured and Derived Properties of Warm H2 in Fields 5 and 6

Region a Filter/BandAssumedFlux Density b AreaBandwidthH2 Line FluxL(0–0 S(3)) nc MH2 d
  Transition(μJy)(arcsec2)(μm)(× 10−18 W m−2)(× 1032 W) (× 106 M)
[1][2][3][4][5][6][7][8][9][10]
F5-R1F1000W0–0 S(3)44 (±7)2.82.02.6 (±0.4)2.8 (±0.4)[4.2] c 5.6 (±1.0)
F5-R2F1000W0–0 S(3)17 (±2)0.62.01.0 (±0.2)1.0 (±0.2)[4.2] c 2.1 (±0.3)
F5-R3F1000W0–0 S(3)7 (±1)0.92.00.4 (±0.1)0.4 (±0.1)[4.2] c 0.9 (±0.2)
F5-R4F1000W0–0 S(3)20 (±3)2.72.01.2 (±0.2)1.3 (±0.2)[4.2] c 2.5 (±0.4)
F5-R5F1000W0–0 S(3)14 (±2)2.62.00.8 (±0.1)0.9 (±0.1)[4.2] c 1.8 (±0.3)
F5-IRSareaF1000W0–0 S(3)87 (±13)30.32.05.2 (±0.8)5.5 (±0.8)4.2 (±0.02)11.0 (±1.9)
F5-IRSLL e 0–0 S(1)30.34.01 e 4.2
F5-IRSLL e 0–0 S(0)30.30.50 e 0.5
F6-AllF1000W0–0 S(3)58 (±9)11.02.04.3 (±0.6)4.5 (±0.7)[4.5] c 14.4 (±2.5)
F6-R1F1000W0–0 S(3)20 (±3)1.62.01.2 (±0.2)1.2 (±0.2)[4.5] c 4.0 (±0.7)
F6-R2F1000W0–0 S(3)6.0 (±0.9)1.32.00.4 (±0.1)0.4 (±0.1)[4.5] c 1.2 (±0.2)
F6-R3F1000W0–0 S(3)5.4 (±0.8)0.72.00.3 (±0.1)0.3 (±0.1)[4.5] c 1.1 (±0.2)
F6-R4F1000W0–0 S(3)7.3 (±1.1)1.02.00.4 (±0.1)0.5 (±0.1)[4.5] c 1.5 (±0.3)
F6-R5F1000W0–0 S(3)7.7 (±1.2)1.32.00.5 (±0.1)0.5 (±0.1)[4.5] c 1.5 (±0.3)
F6-IRSboxF1000W0–0 S(3)76 (±11)30.32.04.6 (±0.7)4.8 (±0.7)4.50 (±0.02)15.3 (±2.5)
F4-IRSonly f SL-LL0–0 S(0)–0–0 S(5)30.32.12.25.00 (±0.03)14.8 (±2.2)

Notes.

a The regions for F5 (Field 5) and F6 (Field 6) are shown graphically in Figure 12. b The flux density is measured from the surface brightness in the JWST MIRI image over the apertures shown in Figure 12, assuming a pixel size of 0.11 × 0.11 arcsec2. c Power-law index n, assuming the warm H2 column density follows a power law with temperature of the form δ NTn δ T, following the work of Togi & Smith (2016) and as discussed for Stephan's Quintet by Appleton et al. (2017). We assume that each of the regions extracted from the MIRI image has the same power law as the average for the whole region, which is a significant limitation of the current analysis. d Estimated mass in warm H2 for T > 100 K assuming a value for the power-law index, n, for the resolved regions that is the same as that measured by the IRS over a larger 5.5 × 5.5 arcsec2 region. These warm H2 mass estimates for T > 100 K employ the method of Togi & Smith (2016). For Field 6, n was derived from the full mapping with both Long-Low (LL) and Short-Low (SL) IRS modules, described in Appleton et al. (2017). For Field 5, it was measured using Spitzer IRS Long-Low data for 0–0 S(0) and 0–0 S(1), combined with the flux measured for 0–0 S(3) from the MIRI data. Uncertainties in the final masses represent propagation of the formal uncertainties and do not include the unknown effect of possible variations from region to region in the actual power-law index, which we have assumed here to be constant. e Line fluxes measured directly from the extracted IRS spectrum from LL IRS spectrograph shown in Figure 3(b). f Unambiguous H2 emission was not easily extracted from the F1000W data for this region (see text). This entry is taken from the IRS data from Appleton et al. (2017) for comparison with the other regions. Note the large value for the power-law index (n), indicative of cooler H2 temperatures near the northern tip of the SQ filament in the SQ-A region.

Download table as:  ASCIITypeset image

Table 2. Integrated Properties of CO Emission in Fields 5 and 6 over the Same Areas as Table 1

Region Vhelio FWHMΣ(Sv ΔV)Σ(Sv ΔV) a M(H2) a Mgas b
   CO(2–1)CO(1–0)×106 ×106
 (km s−1)(km s−1)(Jy km s−1)(Jy km s−1)(M)(M)
[1][2][3][4][5][6][7]
F5R1A6401 (±6)101 (±15)0.9 (±0.1)0.28 (±0.05)18.8 (±3.0)25.6 (±4.1)
F5R1B6646 (±4)73 (±10)0.6 (±0.1)0.19 (±0.03)13.0 (±2.0)17.7 (±2.7)
F5R26388 (±4)119 (±10)0.8 (±0.1)0.25 (±0.02)16.7 (±1.5)22.7 (±2.0)
F5R36656 (±2)105 (±5)1.5 (±0.1)0.47 (±0.02)31.5 (±1.5)42.8 (±2.0)
F5R46618 (±4)84 (±8)1.1 (±0.1)0.33 (±0.04)22.4 (±2.4)30.4 (±3.3)
F5R56537 (±5)67 (±12)0.6 (±0.1)0.20 (±0.04)13.5 (±2.7)18.3 (±3.6)
F6-ALL6054 (±2)91 (±4)4.7 (±0.2)1.5 (±0.1)99.1 (±4.4)134.8 (±6.0)
F6R16073 (±1)47 (±2)1.1 (±0.1)0.3 (±0.0)23.0 (±1.1)31.3 (±1.4)
F6R26040 (±2)93 (±4)1.2 (±0.1)0.4 (±0.0)26.1 (±1.3)35.4 (±1.7)
F6R36051 (±2)41 (±5)0.2 (±0.0)0.1 (±0.0)4.3 (±0.6)5.9 (±0.8)
F6R46009 (±2)37 (±4)0.3 (±0.0)0.1 (±0.0)5.9 (±0.7)8.1 (±1.0)
F6R56103 (±4)48 (±10)0.2 (±0.1)0.1 (±0.0)4.9 (±1.1)6.7 (±1.5)

Notes.

a We assume a conversion between a line flux at CO (1–0) to that of CO (2–1) of SCO(1−0)/SCO(2−1) = ${({\nu }_{1-0}/{\nu }_{2-1})}^{2}{({r}_{21})}^{-1}$, where ν1−0 and ν2−1 are the rest frequencies of the transitions, and r21 is assumed to be 0.8, similar to that measured in the Taffy Galaxies bridge (Zhu et al. 2007), which shares many similarities with the gas in Stephan's Quintet (see Appleton et al. 2022). The H2 masses presented here are for an XCO value = XCO,20 = N(H2)/ICO = 2 × 1020 cm−2 (K km s−1)−1, the standard value assumed for our Galaxy. In the text we discuss how this may not be applicable in such a turbulent region. MH2 = 7.72 × 103 D2Σ(Sv ΔV)(1+z)−1, where D = 94 Mpc and Σ(Sv ΔV) is the CO(1−0) line integral with Sv in janskys and the velocity V of the gas in kilometers per second. We assume z = 0.02. b Total gas mass Mgas = MH2 × 1.36, includes a 36% correction for helium (Bolatto et al. 2013).

Download table as:  ASCIITypeset image

In Table 1, Column 1 gives the name of the regions defined in Figure 12 for Fields 5 and 6. The regions were chosen to cover the main warm H2 features. Column 2 gives the MIRI filter used or the Spitzer module used (see later). For most of the regions, we extract the equivalent 0–0 S(3) flux from the JWST F1000W image. The pure rotational molecular hydrogen transition is listed in Column 3. Column 4 gives the measured H2 flux density measured from the JWST F1000W over an area in arcsecs given in Column 5. We convert these flux densities to equivalent line fluxes, Column 7, using the filter bandwidth listed in Column 6. The line luminosity of the 0–0 S(3) line is given in Column 8, assuming a distance to SQ of 94 Mpc. We also present the properties of warm H2 taken from Spitzer for the regions, and the appropriate IRS modules and H2 rotational transition is listed. These are used below in the calculation of H2 masses. Uncertainties are provided in parentheses for each measured quantity.

Figure 12.

Figure 12. Regions used to determine the brighter warm H2 flux hot spots and CO spectral extractions areas for (a) Field 6, including five regions (red circles) and the full region (yellow), and (b) Field 5, showing the five extraction regions. Larger regions (not shown) covering an area 5.5 × 5.5 arcsec2 were also extracted, corresponding to the the area sampled by Spitzer IRS. The properties of the extracted fluxes and CO spectral properties are presented in Tables 1 and 2. Contours are of the integrated CO(2–1) emission. Smaller extraction regions concentrating on the properties of the small-scale CO emission for all three regions are presented in Appendix A. The white scale bar is 3'' in length.

Standard image High-resolution image

To calculate the total warm H2 mass in the regions, we make the following assumptions. First, since in most cases we only have the line flux of the 0–0 S(3) line derived from the MIRI image, we must assume something about the temperature distribution of the gas in order to estimate the total warm H2 mass. Since the regions we are interested in were observed with Spitzer over a larger area, but over many of the rotational H2 lines, we can assume, as an approximation, that the temperature distribution measured by the Spitzer IRS (Appleton et al. 2017) can be applied to all the smaller regions resolved by MIRI imaging. In that paper, we showed that the excitation diagrams for many regions extracted over SQ were almost always consistent with a single power-law relationship between the warm H2 column density and the H2 temperature of the form δN(H2) ∝Tn δ T, following the work of Togi & Smith (2016). For Field 6, the Spitzer data provided a power-law index n = 4.5, which is quite shallow compared with normal H ii regions, and consistent with warm gas that is shock heated (Appleton et al. 2017).

Because Field 4 (the region in SQ-A) was ambiguous regarding isolating the H2 line in the F1000W filter because of contamination from neighboring broad wings of PAH features (Figure 3), we did not attempt to measure the masses of smaller regions in this field. Rather we provide, for comparison with Fields 5 and 6, the IRS-measured values. This analysis yielded a much steeper value of n = 5. This is consistent with our observation that Field 4 is dominated by star formation.

Field 5 is a little more complicated. Unlike the other cases, this region lies outside the full spectral mapping region performed by Spitzer, and only the IRS Long-Low module was used to cover this area. This provided measurements of the 0–0 S(0) 28 μm and 0–0 S(1) 17 μm transitions, but not the shorter-wavelength transitions. However, the JWST observations from MIRI can provide the 0–0 S(3) flux, if we extract over the same extraction region (referred to as F5-IRSarea in Table 1) observed by Spitzer. With three points on the excitation diagram given in the table, we were able to measure a power-law index of 4.2. The flat value for n is comparable with some of the warmest regions of Stephan's Quintet (Appleton et al. 2017).

The above measurement assumes that the MIRI and Spitzer data can be compared reasonably accurately given that one is obtained from an image and the other from a spectrometer. To check the validity of this comparison, we compared the measured 0–0 S(3) line flux from IRS with the flux we derived from the extraction of the same area for Field 6. As Table 1 shows, we measured a line flux of 4.6 (±0.7) × 10−18 W m−2 for the 0–0 S(3) line from the JWST image. With the IRS, we measured (region 142 of Appleton et al. 2017) a line flux of 4.7 (±0.3) × 10−18 W m−2. The agreement provides some confidence that we are measuring the fluxes correctly using the MIRI F1000W image.

Having established an average power-law index, we now apply that same value of n appropriate to the smaller extracted regions for Fields 5 and 6 (Column 9 of the table). Without spectroscopic observations on the subarcsec scale, this is obviously an approximation. Using the 0–0 S(3) flux measured by JWST, we estimate the total H2 mass in these regions integrated over the power law down to a temperature of 100 K. These masses (and uncertainties) are given in Column 10 of Table 1. We caution that these uncertainties do not include the unknown effect of possible variations, from region to region, in the actual power-law index, which we have assumed is equal to the average value derived over a much larger area.

From this analysis, we see that the warm H2 mass of the entire arrow-head like structure in Field 5 (F5R1) is 5.6 (±1.0) × 106 M and its core (F6R2) contains almost half of the mass, 2.1 (±0.3) × 106 M. The partial ring structure (F5R3) has a much smaller warm H2 mass of 0.9 (±0.2) ×106 M. The two more distant regions of warm H2 gas at the end of the filament in Field 5 have larger warm gas masses of (F5R4) 2.5 (±0.4) and (F5R5) 1.8 (±0.3) × 106 M, respectively.

In the case of the head-tail structure of Field 6, the triangular hot spot (F6R1) has a warm H2 mass of 4 (± 0.7) × 106 M and most of the clumps in the tail have masses of 1–2 × 106 M. The total H2 mass integrated over the whole structure is 14 (± 3) × 106 M, suggesting that there is considerable extended emission in the tail.

Fields 4, 5, and 6 have similar total warm H2 masses when integrated over the larger 5.5 × 5.5 arcsec2 (2.5 × 2.5 kpc) area measured by Spitzer IRS.

How do these warm H2 masses compare with the mass of molecular gas measured by ALMA? In Table 2, we present the measured and derived properties of the gas from the ALMA observations. We use the integrated CO(2–1) line profile for the same regions as those used in Table 1 for the warm H2. Here, we are concerned with measuring the total CO line profile associated with each of the regions for which we measure warm H2 emission. Column 1 gives the region description. Columns 2 and 3 provide the mean heliocentric velocity and FWHM of the emission after fitting each region with a Gaussian line profile. We note that for F5R1 there we two detected CO emission features, referred to as F5R1A and B. This is because R1 contains emission from both the compact core feature and the linear filament. The inclusion of two distinctly different structures (with different kinematics) within the one region is clearly seen in Figure 12(b). These structures will be discussed in more detail in Appendix A. Column 4 provides the line integral of the CO(2–1) emission in jansky-kilometer per second, and Column 5 provides an estimate of the equivalent CO(1–0) emission, based on an assumed ratio of the integrated intensity of the main-beam temperature in the two lines (r21 = I[12CO(2–1)]/I[12CO(1–0)]). Values of r21 range from 0.76 to 0.9 in normal galaxies at low and intermediate redshift, with a slightly larger range for perturbed galaxies (Braine & Combes 1992; Daddi et al. 2015). We assumed a value of 0.8 to be similar to that measured in the highly turbulent Taffy Bridge (Zhu et al. 2007). Column 6 provides an estimate of the H2 mass for an assumed value of the X-factor, XCO. For the purposes of this discussion, we present the H2 masses in Column 6 and the total gas mass in Column 7, assuming XCO is the average Galactic value. This value is likely to be significantly overestimated in regions of strong shocks and turbulence. The footnotes to the table provide more details about the definitions and assumptions used in the calculations.

Bearing in mind the uncertainties in both the warm and cold molecular masses, we can estimate the ratio of M(H2)warm/ M(H2)cold by comparing the warm gas masses estimated by JWST with those from ALMA. For the F5R1 regions (the arrow-head structure in Field 5) the ratio of the warm-to-cold H2 masses is 5.6/(25.6 + 17.7) = 0.13. This assumes that both the CO linear filament and the CO core both contribute to the total mass. It is not clear that the linear structure is actually detected in warm CO, so the ratio could be much higher if that object was excluded (0.22). However, in reality, since XCO is likely to be much smaller than Galactic (say one-fifth XCO,20), then the fraction of warm gas with T > 100 K is likely to be very large (almost equal to the cold H2). These values are significantly greater than those seen in normal galaxies. F5R3, which includes the partial ring of CO emission, has a warm-to-cold H2 ratio of 0.02, which is much more like a normal galaxy. Regions F5R4 and R5 show warm-to-cold ratios somewhat in between the two extremes (0.08 and 0.1). These regions lie at the opposite end of a warm filament (only weakly detected in CO) that connects them to the other Field 5 structures. The warm triangular hot spot in Field 6 shows an average ratio of warm-to-cold H2 = 0.12. Again, if we assume XCO is one-fifth of the Galactic value, this would imply a significant warm component. F6R2 has a much lower ratio of warm-to-cold gas (0.03), a number which is closer to that of normal star-forming galaxies. It is interesting that this is the region which emits strongly in the F770W filter, which is consistent with this containing a PDR. The other clumps in the tail have unusually large values of warm-to-cold CO (in the range 0.18–0.22), reflecting the relative faintness of the CO in the clumps in the tail of Field 6. The overall conclusion is that, except for the regions which show obvious star formation (PAH emission and compact star clusters embedded in the emission), the warm-to-cold H2 ratios are unusually large for both Field 5 and Field 6. We will discuss this further in the discussion in the context of shocks and turbulence.

We note that our ALMA radio continuum observations at 230 GHz might, in principle, provide useful information about the mass of cool dust associated with the target regions. Inspection of the continuum maps did not detect dust emission from any of the regions of interest down to an rms upper limit of 0.03 mJy beam−1. We assume a dust mass Mdust = Sν D2/(κν B(ν,Tdust)), where Sν is the flux density at frequency ν, D is the distance to the galaxy, and B(ν,Tdust) is the blackbody function at frequency ν and temperature Tdust. Assuming a reasonable average value for the dust absorption coefficient κ1.3 mm of 0.03 m2 kg−1 (extrapolating from 500 μm from Clark et al. 2019), Tdust = 20 K, we estimate a 3σ upper limit to the dust mass of Mdust < 5.6 × 105 M beam−1. A dust temperature of 20 K is consistent with previous Herschel observations of dust emission on much larger scales (Appleton et al. 2013). If the typical gas-to-dust ratio for galaxies (Sandstrom et al. 2013) is 80–100 (this number may be higher if dust is destroyed in the shock waves), the corresponding upper limit to the gas mass is at least a factor of 10 larger than the CO-derived gas masses within the same beam. This suggests that the continuum observations are not deep enough to provide useful physical constraints on the cold dust properties of the observed regions. Observations at higher frequencies would be desirable to potentially detect cold dust from these regions.

5. Discussion

SQ is a fantastic example to study the dynamical interaction between gas phases in a group environment. The big leap forward in spatial resolution offered by the JWST images compared to Spitzer sheds light on the morphological structures of the multiphase gas, and the physical processes at the origin of the formation and excitation of the molecular gas in the IGM. We discuss what we are learning about these two aspects in the following sections.

5.1. Morphology of the Molecular Gas Emission and Large-scale Kinematics

In order to explain the formation of large masses of molecular gas and high luminosity in the pure rotational ground-state transitions of H2, Guillard et al. (2009) presented a model in which a large-scale driving shock wave (traveling at 600–800 km s−1) from the intruder passes into a multiphase medium of gas in a filament of tidal debris caused by previous tidal interactions between group members (see Figure 13 for a schematic view). We summarize here that scenario, which was based on observations with limited spatial coverage and resolution with Spitzer, and then discuss what we are now learning from the new JWST plus ALMA data presented here about the link between the large-scale dynamics in the group and the formation and excitation of the molecular gas in the IGM.

Figure 13.

Figure 13. Sketch of the scenario leading to the formation of multiphase gas in the post-shock region. The relative motion between the intruder, NGC 7318b, and the tidally stripped filament drives a large-scale shock (40 kpc long, symbolized with the red line) in the low-density gas (the light blue filament). The pre-shock medium being multiphase, shock velocities are slower in the denser (atomic) regions (darker blue clouds), allowing molecular gas (dark blue clumps) to form in less than the crushing timescale of the H i clouds (see Section 5.1 and Guillard et al. 2009). We argue in Sections 5.2 and 5.3 that the small-scale (a few 100 pc) H2 and CO structures observed by JWST and ALMA may be the result of the fragmentation of the clouds and subsequent reformation of a "foggy" molecular gas (symbolized as tiny green droplets). At small scales, as sketched in the zoom-in panel on the right, the dynamical interaction between the flow of hot gas (red arrow) and the resulting cloudlets drives a cycling across ISM phases, leading to the formation of extended and turbulent H2 tails.

Standard image High-resolution image

The large-scale shock wave driven by the intruder, NGC 7318b, passes over the tidal, multiphase filament, heating low-density gas (nH < 2 × 10−2 cm−3) to X-ray-emitting temperatures (see O'Sullivan et al. 2009), and destroying any dust grains present there (illustrated in Figure 13 as the light orange medium). Regions in the pre-shock region with overdensities greater than 10 (typically H i gas with nH > 0.1 cm−3) compared to the volume-filling gas are heated to lower (< 106 K) temperatures, and can cool and form molecules on dust grains as they survive the slower shocks (sketched as the dark blue clouds in Figure 13). Dust grains are expected to survive in shocks with velocities lower than 200 km s−1 (Jones & Nuth 2011). Guillard et al. (2010) have indeed shown that dust survives in the denser regions of the IGM of the galaxy collision, and that the PAH and continuum IR emission is consistent with the dust associated with the warm (>150 K) H2 gas, heated by a UV radiation field of intensity comparable to that of the solar neighborhood, G0 = 1.4. Natale et al. (2010) also showed that a diffuse IR dust emission component is associated with the collisional cooling of the hot halo gas, provided that dust is injected into the hot plasma by tidal stripping or stellar feedback. In the next section, we discuss how the dynamical interaction between the cold and hot gas can lead to the same physical conditions downstream from the shock, as thermally unstable gas is produced within the turbulent mixing layers (Gronke & Oh 2022), leading to a cycle of molecular clouds breaking up and reforming in the flow. This is sketched as the turbulent tail in the zoom-in panel on the right-hand side of Figure 13.

An efficient transfer of the bulk kinetic energy of the galaxy collision to turbulent energy within the molecular gas, associated with a continuous cold gas destruction–reformation cycle, is required to make H2 a dominant coolant of the post-shock gas (Guillard et al. 2009). In the next section, we discuss how the spatial decoupling between cold and warm gas at GMC scales, revealed by the comparison between ALMA and JWST images, is shedding light on the physical processes that drive this energy transfer mediated by thermal instability and a dynamical interaction between gas phases.

5.2. Cycling and Redistribution of Matter across Gas Phases: Shattering and Growth of Cold Clouds in the Intergalactic Medium

The ALMA and JWST data presented in this paper reveal complex and turbulent structures, with large velocity dispersions (40–100 km s−1) at GMC scales, and a spatial decoupling of cold and warm molecular gas. Although a quantitative interpretation of these data is difficult and beyond the scope of this paper, we qualitatively discuss a scenario that relate the cold and warm H2 distributions to theoretical work of cold clouds shattering in their dynamical interaction with the hot volume-filling hot plasma and subsequent gas cooling (Gronke & Oh 2020).

When the shock wave hits the multiphase IGM filament, the clouds are compressed on a crushing timescale, tcrush = Lc /Vc , where Vc is the shock velocity in the cloud and Lc the size of the pre-shock cloud. The shocked clouds experience dynamical instabilities and mixing at the interface with the hot gas flow, as well as thermal instabilities induced by gas cooling, which provide an additional fragmentation mechanism occurring over the gas cooling timescale, tcool, as the shock moves into the cloud. The gas stripped from the cloud may cool and condense as a result of thermal instability. This occurs if the gas cooling time tcool is smaller than tcrush. There is a wide range of parameter space (typically nH ≈ 0.1–100 cm−3, Lc = 1−100 pc) where tcrush > tcool, i.e., molecular gas can form before the shock has crossed the cloud and thus before dynamical instabilities fully develop (see Figure 4 in Guillard et al. 2009). The cooling time in the warm molecular cloudlets is mostly driven by H2 cooling, therefore the dust-to-gas mass ratio could, in principle, indirectly affect the H2 reformation in the post-shock gas, and therefore the cooling time. However, modeling of the dust-to-gas ratio, and the the PAHs to very small grains abundance ratio, based on Spitzer and Herschel spectroscopy and photometry, show that these are comparable to the Galactic value in the molecular gas (Guillard et al. 2010). This again suggests that dust survives the slower shocks driven by the dynamical instabilities and mixing into the molecular cloudlets.

Because of thermal instability, the post-shock clouds shatter in many clumplets, which then are seeds for further cooling-induced condensation, but also mixing (Farber & Gronke 2022a). In this parameter space where the cooling time is shorter than the cloud destruction timescale, this shattering process causes the gas to cycle through the ISM phases with the reformation of cold clumpy gas. Cloud growth depends on density as well as on local heating processes (Jennings et al. 2023). Given the large masses of diffuse molecular gas observed in the IGM where atomic H i gas is undetected, qualitatively, the post-shock physical conditions must be in the density and irradiation regime where the cloud evolution is dominated by shattering and perhaps cooling-driven coagulation.

Many of the structures detected with JWST along the main north–south molecular filament are shaped like cometary tails, with Field 6 being a prime example (see Figure 12). They could be signatures of a misty H2 outflow driven by the dynamical interaction between cold and hot ISM phases, e.g., ram pressure stripping of CO clouds as the hot X-ray-emitting gas flows past them. This is illustrated in the zoom-in panel of Figure 13. The high abundance of warm H2 relative to the CO-derived cold gas mass discussed in Section 4.1 for both Field 6 and Field 5 lends support to this idea, despite their different spatial positions in the SQ group. We will return to these differences later. The pressure-driven flow velocity is larger than the turbulent gas velocity because the droplets are tiny. At the scale of an H2 droplet, the difference between the turbulent velocity dispersion and the flow velocity is even larger. Farber & Gronke (2022a) suggest that in typical circumgalactic medium or galaxy cluster environments, cloud shattering during the cooling process could be responsible for in situ formation of a fog of warm (1000 K) gas embedded in hot (107–8 K) plasma, with molecular shattering length scales of the order of 0.1–100 pc. This would imply that fragments formed in such manner are very small, with a low volume-filling fraction. A detailed investigation of the physical conditions, which could be achieved with JWST H2 spectroscopy, would allow us to determine whether the cloudlets can coagulate or whether the molecular streams we observe remain misty.

5.3. Multiscale Kinematical Structures and Survival of the Cold Gas

The very complex, multiscale kinematical and morphological structures observed at high resolution with JWST and ALMA points toward a highly clumpy (misty) structure of the molecular gas in the IGM of SQ. This could explain why the Lyα line emission is so strong and the line profile so broad (Guillard et al. 2022). Indeed, the cloudlets must have a small volume-filling fraction but their surface filling factor must be close to unity so Lyα photons can escape through multiple scatterings off the clump surfaces without being absorbed in the inter-clump, dust-free plasma.

The ALMA and JWST observations also shed light on a puzzle raised by the H2 observations by Spitzer (Appleton et al. 2006), as well as single-dish CO observations, which showed extreme velocity dispersion (≈1000 km s−1) of the molecular gas on kiloparsec scales. These large-scale velocity gradients, also observed in other gas phases, could partly originate from coherent bulk flows (Guillard et al. 2022). The H2 molecules can only survive at shock velocities ≲20 km s−1, which is much smaller than those velocities at kiloparsec scales, but also smaller than the line widths of the CO clouds. The ALMA observations presented in this paper indeed show that the CO clumps have velocity dispersions in the range of 40–100 km s−1. Those velocities are also higher than the shock velocities needed to account for the excitation of the pure rotational lines of H2 of the order of 5–20 km s−1 (Guillard et al. 2009; Lesaffre et al. 2013; Appleton et al. 2017).

The tiny H2 droplets are unlikely to survive over the length of the tails, which supports the idea that they are continuously destroyed and reformed out of the more diffuse gas. This could explain why some of the CO clumps (in the tail of Field 6, for instance) have high-velocity dispersions (see Figure A1), as they are reformed out of pressure-driven outflow from the parent cloud. This will need to be further investigated with future H2 spectroscopy with JWST to study the velocity gradients of the warm H2 gas at the scales of the ALMA observations. Spectroscopy would also allow us to see whether this fog of small warm clouds is being accelerated in the bulk flow of the hot X-ray-emitting gas, and would likely differ from the motion of the massive CO clouds. The temperature of the warm droplets may also increase along the tail as the foggy clouds progressively mix and become heated in the faster flow. This might be evident in the excitation diagram of H2 when a full complement of H2 rotational lines can be measured. Those observations will motivate a direct comparison with future numerical simulations.

The comparison between the warm and cold molecular gas spatial distributions in the IGM of SQ shows that the powerful mid-IR emission originally discovered with Spitzer spectroscopy does not primarily have a one-to-one correlation with the cold CO(1–0)-emitting gas. This supports the view of the model put forward in Guillard et al. (2009), which involves a turbulent energy cascade and dissipation of mechanical energy within the warm H2 phase, both driven by the dynamical interaction between the hot and cold gas phases. The complex morphology and kinematics of the warm and cold H2 structures shows that this dynamical interaction occurs over a wide range of scales. This is expected from an energy cascade over many orders of magnitude in spatial scales, from tens of kiloparsecs at the injection scale of mechanical energy from the galaxy collision to subparsec dissipation scale through line emission (Guillard et al. 2009).

5.4. Does the Proposed Molecular Cloud Outflow Model Apply across the Whole SQ System?

As we have stressed earlier, the high luminosity in the warm molecular hydrogen lines over a large part of the Stephan's Quintet system requires a common mechanism to convert a fraction of the kinetic energy resulting from the large-scale shock wave into rotational molecular hydrogen via low-velocity shocks (Guillard et al. 2009; Appleton et al. 2017). We have seen that, of the three regions we studied in detail here, the one that most closely fits the proposed cloudy outflow model is the head-tail structure in Field 6, where we seem to be witnessing a possible breakup of a large cold molecular cloud into warmer H2 in both the head and tail.

Although modeling of the actual structure of a cold molecular cloud in a hot wind will require further simulation, it is likely that cometary structures should be common along the main shock wave, and the diversity of structures should be similar to those seen in the cloud-survival sequences of models by Farber & Gronke (2022b). Those model sequences, where clouds are entrained in a fast wind, show that in some cases the cooler head and tail structures can survive for long periods, whereas in other cases, depending on the properties of the clouds, they can dissolve fairly quickly. Indeed, one of the models has a morphology remarkably similar to the shell-like CO and warm H2 structure in the head of Field 6, including the spike feature at the head of the clump. Many of the shapes of the bright H2-emitting clumps seen along the main shock front in Stephan's Quintet (Figure 2) show suggestions of head/tail structures, and have filamentary structure very reminiscent of the modeled shapes seen when the head/tail structures survive. Others have shapes similar to the models where the clouds have lost their heads and are in the process of dissolving. Therefore, it is plausible that our picture of the foggy outflowing cloud picture is applicable to a large number of cloud structures seen in the main shock.

In addition, Figure 2 shows head-tail structures far to the south and west of the main SQ-A region. One of those tails extends over 30''–40'' (10–20 kpc!). These structures may represent the interaction of dense molecular clumps with the outer hot halo of the intruder far downstream of the main shock. More observations of the cold and warm molecular gas and their kinematics will be needed to test this idea.

In the case of our study of the Field 4 region in SQ-A, it is less clear that the cloud outflow model is applicable, although we note that a faint irregular tail is seen extending to the east of the main concentration of gas and stars. Xu et al. (1999) have previously suggested that SQ-A is too young to have formed as a tidal dwarf galaxy, but is more likely a starburst region created as part of the main shocked filament. The H2 gas temperature measured by Spitzer in SQ-A is much lower than other regions in the main filament (Appleton et al. 2017), and the gas may have cooled to the point that stars can form in large quantities there. Although understanding how those stars form is beyond the scope of this paper, we have shown in Guillard et al. (2009) that the large (>50 pc), dense (nH > 100 cm−3) pre-shock clouds are expected to be gravitationally unstable in the hot gas at this pressure (see their Figure 4). Therefore, their evolution will be different than the more diffuse gas, and should trigger star formation. This may explain star formation in SQ-A and other star-forming regions scattered in the IGM. These large pre-existing clouds have crushing times longer than the SQ collision age (believed to be a few tens of millions of years). However, such clouds must be rare, otherwise we would see major bursts of star formation all over the main filament, which is not observed.

The case of Field 5 (in the bridge between NGC 7319 and the main filament) is interesting. As we have shown in Section 4, the morphology of the warm and cold gas is not so obviously a head-tail structure, although there is a warm gas filament running through the region. This region may owe its overall morphology and kinematics to another process—for example, the collision of a dense molecular concentration (the compact core) with a diffuse region near NGC 7319, driving a ring into the target material. The elongation of the compact CO-emitting core along the direction which points toward the center of the CO-emitting ring, and the existence of a clump of CO in the ring at the same velocity as the core, may support a dynamical connection. Future warm H2 spectroscopy will determine whether the warm gas connecting the two CO structures with disparate radial velocities is related, or merely chance alignments of two different cloud systems along the line of sight. Such events, if they occur, are likely to be relatively rare, and are unlikely to provide a generic mechanism for the dissipation of collisional kinetic energy into widespread warm H2 emission.

Finally, we come to the complex network of filaments around NGC 7319, which, for the most part, appear pink in the false-color representation of Figure 2. This implies that the 7.7 μm and 15 μm emission is much stronger here than in the regions which are warm H2 dominated.

The irregular radial structure of the filaments might suggest a nuclear origin connected with the active galactic nucleus (AGN) in NGC 7319. The JWST MRS data in the central region of the AGN indeed reveal a spatially resolved ionized and molecular outflow, with kinematical evidence of an interaction between the AGN jet and the host ISM (Pereira-Santaella et al. 2022). Perhaps previous jet activity from the AGN has somehow formed these strange filaments.

6. Conclusions

We have made a detailed comparison, at subarcsec resolution, between previous HST, recent ERO JWST NIRCam and MIRI imaging, and ALMA CO(2–1) imaging spectroscopy of three targeted regions approximately 3 kpc in size in the large-scale (40–50 kpc) intergalactic gas in Stephan's Quintet. Based on a comparison with a full spectral map by Spitzer and observations made by JWST, we demonstrated that much of the gas along the main molecular filament detected in the F1000W/MIRI filter is dominated by strong emission from the 0–0 S(3) line originating in warm molecular hydrogen. This allowed us to reveal differences in the spatial distribution of the warm H2 gas phase, measured by JWST, and the colder CO-emitting gas, at comparable subarcsec resolution in three regions targeted by ALMA. Our study has come to the following conclusions:

  • 1.  
    The three zoomed-in regions we studied show a great variety of structure in the warm H2 when compared with the cold molecular gas on scales from 150 pc to 3 kpc. In many cases, diffuse ionized gas (observed by HST and F200W/NIRCam) follows the warm H2 distribution in the two structures believed to be shock dominated (ALMA Fields 5 and 6), whereas the cold molecular gas was more clumpy and not always strongly correlated with the warm H2. Both structures show an overabundance of warm H2 over the colder gas with reasonable assumption about the conversion from CO line intensity to total H2 mass. In the star-formation-dominated region (at the center of ALMA Field 4), the 10 μm and 7.7 μm bands are well correlated, as expected in regions dominated by PDRs.
  • 2.  
    In the (ALMA Field 6) region observed near the center of the main shocked filament, the warm H2 emission resembles a head-tail structure. Strong CO emission forms a partial shell of cold H2, which is embedded in a bright extended hot spot of warm H2 (F1000W emission), containing approximately 4 × 106 M of warm gas if we assume the excitation of the gas is similar to that inferred on a larger scale by Spitzer. This triangular-shaped head also shows extended ionized gas. The tail is composed of clumps of warm and cold H2 emission scattered in two streams eastwards from the head. Near the apex of the triangular warm H2 head and in some regions along the tail, CO line widths of 60 < FWHM (km s−1) < 100 km s−1 are observed, suggesting quite turbulent (probably gravitationally unbound) gas. We suggest that we are witnessing the shattering of a dense cold molecular gas cloud in a fast hot wind, leading a fog of warm H2 (see below). We believe that this process is a commonly occurring phenomenon in Stephan's Quintet, judging by the existence of other head-tail strcutures along the main shock and elsewhere in the group.
  • 3.  
    In another shock-dominated warm H2 region (ALMA Field 5), which is not part of the main molecular filament, we notice a different morphology. Warm molecular gas appears to connect together a compact elongated CO-emitting clump that points toward the center of a ring of CO emission separated by 1farcs5 (∼700 pc). The two CO structures have very different radial velocities, with the compact core having strongly blueshifted emission with respect to the ring. Both the ring and the core show broad (∼100 km s−1 FWHM) CO line widths, suggesting that they are highly turbulent. One possibility is that a dense molecular cloud (perhaps from the outer parts of the intruder galaxy) has impacted a low-density gas filament in the outer parts of NGC 7319, entraining gas and dissipating kinetic energy, which heats the bridge between them. Such events are probably rare. Alternatively, the connection between the two CO clouds by the warm gas may be entirely coincidental.
  • 4.  
    The distribution of warm, cold, and ionized gas in the SQ-A star-forming region (ALMA Field 4) is more typical of the gas distribution of a small dwarf galaxy. The warm and cold molecular gas are contained within a roughly elongated stellar distribution (as observed by NIRCam). The lack of clear rotation in the CO-emitting gas may argue against this object being created in a tidal stream. Instead, we suggest it may be more consistent with the collapse of a rare, pre-existing dense molecular structure which has become gravitationally unstable after it was compressed by the passage of the large-scale shock, leading to a starburst.

We present a conceptual picture for the transfer of kinetic energy from the intruder galaxy's large-scale shock into a pre-existing multiphase tidal filament via the formation and subsequent molecular cloud shattering of cold molecular clouds in a post-shock layer. After the main shock passes over a multiphase tidal filament, both X-ray-emitting hot gas and dense molecular clouds can form together if the molecular formation time on dust grains is shorter than the cloud-crushing time of the clouds. However, as the molecular clouds grow, the velocity of the denser forming clouds relative to the diffuse hot gas causes them to experience ram pressure stripping from the hot component. Models show that, in this situation, the cold molecular clouds can shatter into a fog of tiny molecular clouds, heated by an exchange in energy with the hot medium. The fog clouds are expected to have low internal velocity dispersion, which we associate with the clouds that radiate mid-IR H2 emission (low shock velocities in the molecular gas). The warm and cold gas associated with ALMA Field 6 is likely one example of many such cold clouds being turned into a warm H2 fog all along the main SQ large-scale shock front. Regions outside the main molecular shock front also seem to show a head-tail morphology (for example, one extremely long, 20 kpc tail to the southwest of the SQ-A star-forming region). These may be examples of a similar process, where the flow of hot gas may originate in the relative motion of the intruder galaxy's hot halo relative to possible pre-existing molecular clouds in the system. Future mid-IR spectroscopy will allow us to test this picture, as the kinematics of the warm H2-emitting fog is likely to be very different from that of the cold clouds, as they rapidly accelerate up to the flow velocity of the hot gas before eventually being destroyed.

Acknowledgments

We acknowledge the remarkable dedication of the scientists and engineers who made the JWST possible, and especially the Early Release Observation (ERO) science team, upon which some of our paper is based. This work is based, in part, on observations made with the NASA/ESA/CSA JWST. The data were obtained from the Mikulski Archive for Space Telescopes at the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS 5-03127 for JWST. These observations are associated with program #2732. This paper makes use of ALMA data from program ID: 2015.1.000241.S (PI: P. Guillard). ALMA is a partnership of ESO (representing its member states), the NSF (USA), and NINS (Japan), together with NRC (Canada), MOST and ASIAA (Taiwan), and KASI (Republic of Korea), in cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by the ESO, AUI/NRAO, and NAOJ. The National Radio Astronomy Observatory is a facility of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc. E.O'S. acknowledges support from NASA through Chandra Award Number GO8-19112A. P.G. would like to thank the Sorbonne University, the Institut Universitaire de France, the Centre National d'Etudes Spatiales (CNES), the "Programme National de Cosmologie and Galaxies" (PNCG) and the "Physique Chimie du Milieu Interstellaire" (PCMI) programs of CNRS/INSU, with INC/INP cofunded by CEA and CNES, for their financial support. We would like to thank an anonymous referee for very helpful comments on the manuscript.

Appendix A: ALMA and CARMA CO Kinematics of the ALMA Field Regions 4, 5, and 6

In this appendix, we present Tables A1, A2, and A3, which provide observed and derived properties of individually extracted spectra for Field 6, Field 5, and Field 4, respectively, based on the CO(2–1) ALMA data. The regions (e.g., F6reg1) refer to the defined region extracted from each of the ALMA fields shown in Figures A1, A2, and A3, respectively. The figures show a gray-scale representation of the surface density map of the CO integrated emission with elliptical (blue) beam-sized extraction areas marked. The spectra and fitted Gaussian profiles are also shown. At the bottom of Figures A1 and A2, we present the integrated emission profiles. For Figure A3, the integrated spectrum is F4p1.

Figure A1.

Figure A1. Field 6 extracted spectra for CO(2–1) spectra covering the CO-beam-sized extraction regions (blue ellipses on the inset) labeled in the inset at the center of the figure. The inset shows a gray-scale representation of the CO integrated emission. The cyan polygons around the emission show the extracted areas which were summed for an integrated spectrum over the whole region. Gaussian fits to the CO line profiles are shown superimposed on the spectra with a red line. Region 21 shows faint emission and was not fit. The red and blue vertical lines indicate the average velocity of the intruder galaxy NGC 7318b and the average group velocity of the main SQ members, respectively. The integrated spectrum obtained over the cyan polygons is also presented as the last spectrum. The integrated properties of all the spectra are given in Table A1.

Standard image High-resolution image
Figure A2.

Figure A2. Field 5 extracted spectra for CO(2–1) spectra covering the CO-beam-sized extraction regions (green ellipses) labeled in the inset at the center of the figure. The inset shows a gray-scale representation of the CO integrated emission, with contours of the same emission superimposed in red. Gaussian fits to the CO line profiles are shown overlaid on each of the spectra with a red line. Regions 11 and 12 are too faint to fit. The red and blue vertical lines indicate the average velocity of the intruder galaxy NGC 7318b and the average group velocity of the main SQ members, respectively. The integrated spectrum obtained over the cyan polygons is also presented as the last spectrum. The integrated properties of all the spectra are given in Table A2. Note that Regions 1, 2, 3, and 4 have significantly different lower velocities than the rest of the emission. Region 14 also shows a double profile (labeled A and B), with faint emission from the low-velocity structure as well as the higher-velocity structure (13, 14, 15, and 16).

Standard image High-resolution image
Figure A3.

Figure A3. Field 4 extracted spectra for CO(2–1) spectra covering the CO extraction regions (green ellipses) labeled in the inset at the center of the figure. The inset shows a gray-scale representation of the CO integrated emission, with contours of the same emission superimposed in red. Gaussian fits to the CO line profiles are shown overlaid on each of the spectra with a red line. The red and black vertical lines indicate the average velocity of the intruder galaxy NGC 7318b and the average group velocity of the main SQ members, respectively.

Standard image High-resolution image

Table A1. Properties of CO Emission for Beam-sized Extractions in Field 6 (see Figure A1)

Region Vhelio FWHMΣ(Sv ΔV)Σ(Sv ΔV) a M(H2) a Mgas b
   CO(2–1)CO(1–0)× 106 × 106
 (km s−1)(km s−1)(mJy km s−1)(mJy km s−1)(M)(M)
F6reg16052 (±1)36 (±3)92 (±8)29 (±2)1.9 (±0.2)2.6 (±0.2)
F6reg26067 (±2)27 (±5)38 (±7)12 (±2)0.8 (±0.1)1.1 (±0.2)
F6reg36069 (±2)40 (±4)98 (±10)31 (±3)2.0 (±0.2)2.8 (±0.3)
F6reg46071 (±1)46 (±3)113 (±9)35 (±3)2.4 (±0.2)3.2 (±0.3)
F6reg56065 (±2)56 (±6)88 (±10)28 (±3)1.8 (±0.2)2.5 (±0.3)
F6reg66082 (±2)50 (±4)113 (±10)35 (±3)2.4 (±0.2)3.2 (±0.3)
F6reg76086 (±6)102 (±13)85 (±12)27 (±4)1.8 (±0.2)2.4 (±0.3)
F6reg86068 (±2)34 (±4)78 (±9)24 (±3)1.6 (±0.2)2.2 (±0.3)
F6reg96072 (±2)42 (±4)75 (±8)23 (±3)1.6 (±0.2)2.1 (±0.2)
F6reg106068 (±2)37 (±4)68 (±8)21 (±3)1.4 (±0.2)1.9 (±0.2)
F6reg116058 (±1)43 (±3)99 (±8)31 (±3)2.1 (±0.2)2.8 (±0.2)
F6reg126061 (±3)77 (±7)96 (±10)30 (±3)2.0 (±0.2)2.7 (±0.3)
F6reg136051 (±3)66 (±6)107 (±11)34 (±4)2.2 (±0.2)3.1 (±0.3)
F6reg146053 (±1)41 (±3)113 (±9)35 (±3)2.4 (±0.2)3.2 (±0.3)
F6reg156021 (±3)91 (±7)156 (±14)49 (±4)3.3 (±0.3)4.4 (±0.4)
F6reg166013 (±2)79 (±6)144 (±11)45 (±4)3.0 (±0.2)4.1 (±0.3)
F6reg176014 (±3)52 (±7)64 (±10)20 (±3)1.3 (±0.2)1.8 (±0.3)
F6reg186007 (±1)20 (±3)42 (±6)13 (±2)0.9 (±0.1)1.2 (±0.2)
F6reg196007 (±1)30 (±3)72 (±8)23 (±3)1.5 (±0.2)2.1 (±0.2)
F6reg206012 (±3)56 (±6)88 (±11)27 (±3)1.8 (±0.2)2.5 (±0.3)
Integrated c 2380 (±230)743 (±70)49.69 (±5.0)67.6 (±7.0)

Notes.

a We assume a conversion between a line flux at CO (1–0) to that of CO (2–1) of SCO(1−0)/SCO(2−1) = ${({\nu }_{1-0}/{\nu }_{2-1})}^{2}$ ${({r}_{21})}^{-1}$, where ν1−0 and ν2−1 are the rest frequencies of the transitions, and r21 is assumed to be 0.8, similar to that measured in the Taffy Galaxies bridge (Zhu et al. 2007), which shares many similarities with the gas in Stephan's Quintet (see Appleton et al. 2022). The H2 masses presented here are for an XCO value = XCO,20 = N(H2)/ICO = 2 × 1020 cm−2 (K km s−1)−1, the standard value assumed for our Galaxy. In the text, we discuss how this may not be applicable in such a turbulent region. MH2 = 7.72 × 103 D2Σ(Sv ΔV)(1+z)−1, where D = 94 Mpc and Σ(Sv ΔV) is the CO(1−0) line integral with Sv in janskys and the velocity V of the gas in kilometers per second. We assume z = 0.02. b Total gas mass Mgas includes a 36% correction for helium (Bolatto et al. 2013). c Integral properties of the whole structure are based on the summed extracted spectrum obtained from the cyan polygons shown in the inset to Figure A1 and the last spectrum in the same figure.

Download table as:  ASCIITypeset image

Table A2. Properties of CO Emission for Beam-sized Extractions in Field 5 (see Figure A2)

Region Vhelio FWHMΣ( Sv ΔV)Σ( Sv ΔV) a M(H2) a Mgas b
   CO(2–1)CO(1–0)× 106 × 106
 (km s−1)(km s−1)(mJy km s−1)(mJy km s−1)(M)(M)
F5reg16416 (±10)145 (±23)0.10 (±0.02)0.030 (±0.005)2.0 (±0.3)2.7 (±0.5)
F5reg26399 (±4)117 (±10)0.15 (±0.01)0.045 (±0.004)3.0 (±0.3)4.1 (±0.4)
F5reg36388 (±5)108 (±12)0.14 (±0.02)0.045 (±0.005)3.0 (±0.4)4.1 (±0.5)
F5reg46412 (±6)75 (±14)0.06 (±0.01)0.018 (±0.004)1.2 (±0.2)1.7 (±0.3)
F5reg56631 (±4)102 (±10)0.15 (±0.02)0.047 (±0.005)3.2 (±0.3)4.3 (±0.4)
F5reg66651 (±4)101 (±8)0.14 (±0.01)0.045 (±0.004)3.0 (±0.3)4.1 (±0.4)
F5reg76640 (±3)107 (±8)0.17 (±0.01)0.052 (±0.004)3.5 (±0.3)4.7 (±0.4)
F5reg86666 (±2)60 (±6)0.10 (±0.01)0.030 (±0.003)2.0 (±0.2)2.7 (±0.3)
F5reg96679 (±2)61 (±6)0.09 (±0.01)0.029 (±0.003)1.9 (±0.2)2.6 (±0.3)
F5reg106688 (±12)177 (±29)0.12 (±0.04)0.037 (±0.014)2.5 (±0.8)3.4 (±1.2)
F5reg11 d
F5reg12 d
F5reg136628 (±6)79 (±15)0.07 (±0.01)0.021 (±0.004)1.4 (±0.3)1.9 (±0.4)
F5reg14A e 6419 (±9)87 (±21)0.052 (±0.02)0.016 (±0.008)1.1 (±0.6)1.5 (±0.8)
F5reg14B e 6671 (±5)59 (±12)0.05 (±0.01)0.015 (±0.003)1.0 (±0.2)1.4 (±0.3)
F5reg156671 (±5)59 (±12)0.05 (±0.01)0.015 (±0.003)1.0 (±0.2)1.4 (±0.3)
F5reg166680 (±3)48 (±6)0.07 (±0.01)0.021 (±0.003)1.4 (±0.2)1.9 (±0.3)
Integrated2300 (±230)718 (±70)48.1 (±4)6.54 (±7)

Notes.

a We assume a conversion between a line flux at CO (1–0) to that of CO (2–1) of S CO(1−0) /S CO(2−1) = ${({\nu }_{1-0}/{\nu }_{2-1})}^{2}$ ${({r}_{21})}^{-1}$, where ν1−0 and ν2−1 are the rest frequencies of the transitions, and r21 is assumed to be 0.8, similar to that measured in the Taffy Galaxies bridge (Zhu et al. 2007), which shares many similarities with the gas in Stephan's Quintet (see Appleton et al. 2022). The H2 masses presented here are for an XCO value = XCO,20 = N(H2)/ICO = 2 × 1020 cm−2 (K − km s−1)−1, the standard value assumed for our Galaxy. In the text, we discuss how this may not be applicable in such a turbulent region. MH2 = 7.72 × 103 D2Σ(Sv ΔV)(1+z)−1, where D = 94 Mpc and Σ(Sv ΔV) is the CO(1−0) line integral with Sv in janskys and the velocity V of the gas in kilometers per second. We assume z = 0.02. b Total gas mass includes a 36% correction for helium (Bolatto et al. 2013). c Integral properties of the whole structure are based on the summed extracted spectrum obtained from the cyan polygons shown in the inset to Figure A2 and the last spectrum in the same figure. d Low signal-to-noise region. e Two features present: 14A = faint broad emission, 14B = brighter higher-velocity emission.

Download table as:  ASCIITypeset image

Table A3. Properties of CO Emission for Beam-sized Extractions in Field 4 (see Figure A3)

Region Vhelio FWHMΣ( Sv ΔV)Σ( Sv ΔV) a M(H2) a Mgas b
   CO(2–1)CO(1–0)×106 ×106
 (km s−1)(km s−1)(mJy km s−1)(mJy km s−1)(M)(M)
F4reg26708 (±1)40 (±3)0.53 (±0.04)0.166 (±0.013)11.1 (±0.9)15.1 (±1.2)
F4reg36731 (±1)37 (±2)0.59 (±0.04)0.183 (±0.012)12.3 (±0.8)16.7 (±1.1)
F4reg46706 (±1)43 (±3)0.19 (±0.02)0.058 (±0.005)3.9 (±0.3)5.3 (±0.5)
F4reg56702 (±1)35 (±3)0.15 (±0.02)0.048 (±0.005)3.2 (±0.3)4.4 (±0.5)
F4reg66731 (±1)35 (±3)0.17 (±0.02)0.054 (±0.006)3.6 (±0.4)5.0 (±0.5)
F4reg76736 (±1)35 (±2)0.21 (±0.02)0.066 (±0.005)4.4 (±0.3)6.0 (±0.4)
F4reg1(int) c 6719 (±1)47 (±1)1.4 (±0.04)0.422 (±0.012)28.2 (±0.8)38.4 (±1.1)

Notes.

a We assume a conversion between a line flux at CO (1–0) to that of CO (2–1) of SCO(1−0)/SCO(2−1) = ${({\nu }_{1-0}/{\nu }_{2-1})}^{2}$ ${({r}_{21})}^{-1}$, where ν1−0 and ν2−1 are the ratios of the rest frequencies of the transitions, and r21 is assumed to be 0.8, similar to that measured in the Taffy Galaxies bridge (Zhu et al. 2007), which shares many similarities with the gas in Stephan's Quintet (see Appleton et al. 2022). The H2 masses presented here are for an XCO value = XCO,20 = N(H2)/ICO = 2 × 1020 cm−2 (K − km s−1)−1, the standard value assumed for our Galaxy. In the text, we discuss how this may not be applicable in such a turbulent region. MH2 = 7.72 × 103 D2Σ(Sv ΔV)(1+z)−1, where D = 94 Mpc and Σ(Sv ΔV) is the CO(1−0) line integral with Sv in janskys and the velocity V of the gas in kilometers per second. We assume z = 0.02. b Total gas mass Mgas = MH2 × 1.36 includes a 36% correction for helium (Bolatto et al. 2013). c Integral properties of the whole structure are based on the summed extracted spectrum obtained from position 1 (see Figure A3).

Download table as:  ASCIITypeset image

In each table, we provide the following information. Column 1 is the region name, Column 2 is the heliocentric velocity of the emission, with uncertainty, Column 3 is the FWHM (and uncertainty) of the spectrum in kilometers per second, and Columns 4 and 5 provide the CO line integral in janskys per kilometer per second for the observed CO(2–1) emission profile and the derived equivalent CO(1–0) line integral, respectively. The assumption for the conversion to a CO(1–0) line integral is given in the footnotes. Column 6 gives the total (cold) H2 mass in each extracted region, assuming a Galactic value for the conversion from CO line intensity to molecular column density (see footnotes). Column 7 gives the total gas mass assuming 36% helium content.

In the shock-dominated regions Fields 5 and 6, the typical mass within the ALMA beam within the brighter regions is a few × 106 M, whereas the brightest regions in Field 4 (SQ-A) are a factor of 8 to 10 times higher, indicating significantly higher hydrogen columns in SQ-A. This difference may be even higher if the value for XCO is larger in the shock-dominated regions.

Appendix B: Photometry and Age Determination from the Fedotov (2014) Catalog

Photometry of the Stephan's Quintet system was performed by Fedotov (2014) using UBVm VI photometry based on observations made using HST and the F336W (U), F438W (B), F547W (Vm ), F606W (V), and F814W (I) filters with WFC3/UV. The magnitudes were converted to the Vega system assuming zero points of 23.484 (U), 24.974 (B), 24.748 (Vm ), 25.987 (V), and 24.680 (I). Corrections for foreground extinction assumed A336 = 0.353 mag, A438 = 0.288 mag, A606 = 0.197 mag, and A814 = 0.122 mag. SSCs were extracted from the images using DAOfind. At the distance of Stephan's Quintet, the WFC3 resolution means that objects as large as 12 pc would be classified as SSCs, but may include more than one cluster. SSCs were selected with a color cut of B – V < 1.5 or V – I < 1.0 in order to minimize detection of foreground objects. In addition, to be considered a SSC the object had to have a photometric error of <0.3 mag, fullfill a sharpness criterion, and exhibit a goodness-of-fit criterion based on the shape of the object compared with the point-spread function.

Table B1 presents B, Vm , V, I photometry for objects we identify in the U-band images, which correspond to the youngest objects associated with the three studied fields. For Field 4 (SQ-A) most of the objects with U-band detections were of moderate age (F4A–F) ranging from 72 to 360 Myr. NIRCam detects a large number of potential star clusters in the body of this likely forming dwarf galaxy, for which we defer discussion to another paper. In Field 5, only one object (F5D) is clearly detected in the U band, with an age of 8 Myr. This SSC is one of the objects that lies in the string of potentially background galaxies to the east and north of the CO ring. In Field 6, F6A and B are the bright young stellar objects (3–5 Myr) with young ages seen associated with the PAH region in the northwest CO clump structure seen in Figure 5(c). The other clumps (F6C, F6D, and F6E) are quite close to the triangular head of the warm H2 emission, and also exhibit young ages. The R.A. and decl. positions are from the original catalog, and were found to be shifted slightly (0farcs25) relative to the Gaia coordinate system used in the paper.

Table B1.  UBVmI Photometric Properties of Clumps in Fields 4, 5, and 6 from the Catalog by Fedotov (2014) a

NameR.A.Decl. UF336W BF435W VmF547M IF814W E(BV) Av Age
 (J2000)(J2000)(mag (unc))(mag (unc))(mag (unc))(mag (unc))(mag)(mag)×107 yr
 [1][2][3][4][5][6][7][8][9]
F4A338.9946933.9804024.371 (0.035)24.758 (0.05)24.253 (0.042)23.1780 (0.032)0.421.36.5
F4B338.9949033.98026724.543 (0.085)24.732 (0.08)23.961 (0.075)22.864 (0.091)0.501.5511.0
F4C338.9949633.98037722.887 (0.066)23.258 (0.07)22.743 (0.076)21.7290 (0.111)0.401.247.2
F4D338.9953933.98025925.552 (0.046)25.234 (0.05)24.322 (0.041)22.805 (0.031)0.501.5536.1
F4E338.9953633.9806924.911 (0.041)24.902 (0.05)24.119 (0.032)22.686 (0.069)0.501.5516.1
F4F338.9956433.98052622.891 (0.14)23.523 (0.12)22.675 (0.116)22.588 (0.149)36.2
F4G338.9960633.98062522.373 (0.03)22.982 (0.05)22.378 (0.055)21.534 (0.039)0.341.057.20
F5D339.0090033.97407527.17 (0.12)28.58 (0.278)28.32 (0.27)26.70 (0.23)0.280.860.79
F6A338.9991833.97287424.70 (0.08)26.04 (0.09)26.10 (0.08)25.91 (0.09)0.080.240.53
F6B338.9991533.97285824.41 (0.10)25.62 (0.08)25.74 (0.10)25.91 (0.13)0.100.310.35
F6C338.9989633.97274826.86 (0.14)28.05 (0.20)27.62 (0.11)27.065 (0.16)0.51.550.3
F6D338.9989633.97258827.028 (0.15)28.31 (0.17)27.914 (0.17)27.711 (0.25)0.51.550.28
F6E338.9991533.97257227.225 (0.17)28.33 (0.19)28.594 (0.28)28.395 (0.39)0.060.180.52

a https://ir.lib.uwo.ca/etd/2635/

Download table as:  ASCIITypeset image

The method used by Fedotov (2014) to estimate the intrinsic Av and age of the clusters was by comparison of the spectral energy distribution (SED) with population-synthesis models. The use of the U band and the V574 filter helps to break the age–reddening degeneracy, as well as avoiding significant contamination with Hα emission in the F606W filter. The age and extinction of each cluster was found by least-squares fitting the SEDs to predictions from (unpublished) updates to the Bruzual & Charlot (2003) models (known as CB07), and made available on request by the authors (S. Charlot and G. Bruzual). These models incorporate the work of Marigo & Girardi (2007) into the evolutionary code, and are described in more detail by Bruzual (2011). They assume a Chabrier (2003) initial mass function. The best-fit values for the E(BV), Av , and age for the clusters are given in the table.

Footnotes

  • 18  

    Other fields were originally put forward in the ALMA proposal, but only these three regions were actually observed.

  • 19  
  • 20  

    We used the IRS 11.3 μm PAH rather than the 7.7 μm PAH, because the signal-to-noise ratio was much better in this PAH feature in the Spitzer data.

  • 21  

    Hereby measured as the FWHM of the emission-line profile, or 2.35× the sigma of these mainly Gaussian lines.

  • 22  

    For example, for clouds of diameter 100 pc (just less than the resolution of ALMA) and given a typical gas surface density of Σgas = 170 M pc−2, the velocity dispersion for a cloud in gravitational (virial) equilibrium would be σ2 = (3/5)π GRcloudΣgas. To be in equilibrium σ = 8.3 km s−1, or a FWHM ∼20 km s−1.

Please wait… references are loading.
10.3847/1538-4357/accc2a